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NOTATION 

Integration  surface  equivalent  to  foil  surface 
2 

Aspect  ratio;  b  /S 

Refers  to  an  airfoil  or  hydrofoil  mean  line  producing  a 
uniform  load  distribution  forward  of  a  and  a  distribution 
decreasing  linearly  to  zero  aft  of  a;  a  is  given  in  fraction 
of  chord 

Coefficient  of  series  representing  pressure  jump  Ap 


Span  of  submerged  foil  section,  measured  perpendicular  to 
free  stream 

Fractional  distance  along  submerged  foil  section,  measured 
perpendicular  to  free  stream 


C,  ,  C  n  p 
La  £6 


CV  C1 


Lift  coefficient;  increases  with  increasing  angle  of 
attack;  lift/qS 

Rate  of  change  of  or  with  respect  to  angle  of  attack 

a  or  6,  respectively,  in  units  of  radians  ^  unless  other¬ 
wise  noted 

Local  lift  coefficient;  (lift  per  unit  span)/qc 


Pitching  moment  coefficient;  positive  tending  to  increase 
angle  of  attack;  (pitching  moment) /qSc' 


CMa’  CM6 


Rate  of  change  of  C„  with  respect  to  angle  of  attack  a  or  6, 
M  -1 

respectively,  in  units  of  radians 


Local  pitching  moment  coefficient;  (pitching  moment  per 
unit  span)/qSc' 

Roll  moment  coefficient;  (roll  moment)/qSc' 


Chord  length  of  lifting  surface,  measured  in  streamwise 
direction 

Mean  geometric  chord;  S/b 

Root  chord  of  section  (chord  at  b'  =  0) 


Location  of  foil  surface 


♦*  Vi 

i  spanwise  pressure  mode 

Depth  beneath  undisturbed  free  surface 
Number  of  chordwise  control  points 
Kernel  function 

Reduced  frequency,  u)c'/2U;  also,  number  of  spanwise  control 
points 

Mach  number 

2 

Stagnation  pressure;  pU  /2 
Submerged  foil  area 
Time 

Free  stream  velocity 

Magnitude  of  orbital  velocity  due  to  waves 

Downwash,  or  perturbation  velocity  induced  normal  to  the 
undisturbed  flow  by  the  foil,  positive  in  direction  opposite 
to  direction  of  positive  lift 

Frame  of  reference  coordinate  parallel  to  flow  and  positive 
in  the  same  direction  as  the  free  stream  velocity 

Unsteady  center  of  pressure;  ratio  of  magnitude  of  unsteady 
moment  in  phase  with  lift  to  the  unsteady  lift  magnitude, 
expressed  as  percent  of  chord  aft  of  the  foil  leading  edge 

Spanwise  position  or  frame  of  reference  coordinate 

Frame  of  reference  coordinate 

Angle  of  attack  relative  to  flow 

Angular  position  in  flapping  motion 

Unsteady  linear  displacement  of  lifting  surface  from  its 
mean  position,  normal  to  surface  at  zero  angle  of  attack; 
also,  angle  of  attack  of  hydrofoil  relative  to  pod 

Dummy  integration  coordinate  used  for  integration  over  z 

Wave  height  from  trough  to  crest 

Dummy  integration  coordinate  used  for  integration  over  y 


xii 


0 


A 

X 

5 

p 

T 

<t> 

w 

0)' 


Dummy  integration  coordinate  used  for  chordwise  integration 

Sweep  angle  of  quarter-chord  line 

Wavelength 

Dummy  integration  coordinate  used  for  integration  over  x 
Fluid  mass  density 

Taper  ratio;  ratio  of  tip  chord  to  root  chord 

Phase  angle  of  loading  relative  to  displacement 

Circular  frequency  of  oscillation  in  rad/sec 

Circular  frequency  of  wave  in  fixed  coordinate  system  in 
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ABSTRACT 

To  provide  guidance  for  hydrofoil  designers,  an 
extensive  operational  evaluation  was  made  of  a  computer- 
based  lifting-surface  theory  for  calculating  hydrodynamic 
loading  on  hydrofoils  in  steady  and  unsteady  motion  in 
inviscid,  subcavitating  flow  at  finite  depth.  Both  the 
numerical  stability  and  the  accuracy  in  comparison  with 
model  data  were  estimated.  By  use  of  an  empirical  cor¬ 
rection,  lift  predictions  can  be  made  with  an  accuracy 
of  about  15  to  20  percent,  while  pitching  moment  pre¬ 
dictions  are  substantially  less  accurate.  Determination 
of  the  accuracy  was  hindered  by  insufficient  or  impre¬ 
cise  data.  Documentation  for  the  computer  program  is 
given. 

ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  partially  funded  under  the  Hydrofoil 
Development  Program  of  the  Naval  Sea  Systems  Command,  Subproject  S4606, 

Task  1703,  and  performed  under  Work  Unit  1153-003.  The  work  was  subse¬ 
quently  expanded  under  the  authority  of  the  Naval  Material  Command  (NAVMAT 
087),  Program  Element  62543N,  Task  Area  ZF-43-421-01,  and  administered  by 
the  Ship  Performance  Department  High  Performance  Vehicles  Program  (1507) 
under  Work  Unit  1500-102. 

INTRODUCTION 

This  report  is  an  evaluation  of  one  version  of  a  commonly-used  general 

approach  (lifting-surface  theory)  to  calculating  fluid-dynamic  loading  on 

airfoils  and  hydrofoils.  This  particular  treatment  was  developed  by  Ashley 
1*  2 

et  al.  and  Widnall  and  programmed  for  computer  use  by  Widnall.  The 

program  will  be  referred  to  as  the  Widnall  program  after  its  author.  The 

theory  is  capable  of  calculating  both  steady  and  unsteady  loading  on  non- 

planar  hydrofoil  systems  at  finite  depth.  It  extends  the  theory  developed 
3 

by  Watkins  et  al.  to  include  nonplanar  surfaces,  finite  depth  effects, 
and  cavitation. 

The  Navy's  interest  in  this  type  of  calculation  arises  from  possible 
design  applications  for  hydrofoil  craft  and  ship  appendages.  These  appli¬ 
cations  include  certain  hydrodynamic  performance  characteristics  (lift  and 


*A  complete  listing  of  references  is  given  on  page  183. 
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side  force,  for  example,  but  not  drag)  and,  in  combination  with  appropriate 
structural  analyses,  predictions  of  hydroelastic  stability  of  lifting- 
surface  structures. 

To  establish  the  usefulness  of  a  calculation  for  design  work,  the  ac¬ 
curacy  of  the  method  must  be  known.  Although  the  Widnall  program  has  been 
available  for  some  time,  a  detailed  analysis  of  its  accuracy  had  not  been 
performed.  Therefore,  the  present  study  was  undertaken.  The  approach  con¬ 
sisted  of  comparing  the  program  predictions  with  a  large  body  of  available 
model  loading  data  from  both  fully-submerged  hydrofoils  and  airfoils.  As  a 
part  of  the  effort,  the  numerical  behavior  of  the  calculation  was  explored 
to  find  the  most  favorable  range  of  numerical  parameter  values  for  program 
operation. 

Some  relatively  limited  evaluations  of  the  theory  have  previously  been 

made.  In  the  initial  development,'*'  a  two-dimensional  version  was  found  to 

4 

agree  with  analytical  results  at  infinite  depth.  Subsequently,  Pattison 
compared  the  two-dimensional  calculation  with  experimental  measurements. 
Significant  differences  were  found,  but  some  of  the  discrepancy  may  have 
been  due  to  the  experimental  configuration  and  measurement  uncertainty. 

In  an  evaluation  of  three-dimensional  loading  predictions,  Langan  and 
Wang’’  compared  results  from  the  Widnall  program  with  experimental  loading 
on  several  airfoils.  Their  work  was  particularly  valuable  in  that  it 
studied  fourteen  other  computer-based  calculations  as  well.  Some  aspects 
of  both  numerical  behavior  and  calculation  accuracy  relative  to  experiment 
were  compared.  None  of  the  fifteen  calculations  was  numerically  stable  to 
three  significant  figures.  Accuracies  were  usually  within  5  percent  of 
airfoil  loading  data,  but  the  Widnall  prediction  of  lift  due  to  camber  was 
somewhat  less  accurate  (within  9  percent).  These  latter  calculations  cor¬ 
responded  to  conditions  of  optimum  numerical  stability  as  determined  by  the 
present  study,  part  of  which  preceded  that  of  Langan  and  Wang. 

These  results  suggest  that  a  program  that  is  more  accurate  than  the 
Widnall  program  can  be  chosen  from  the  above  study.  However,  the  capabili¬ 
ties  of  treating  free-surface  effects,  including  surface-piercing  configu¬ 
rations,  and  unsteady  loading,  not  generally  available  in  the  other 
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programs,  make  the  Widnall  program  potentially  more  useful  in  designing 
hydrofoils.  Therefore,  the  present  study  is  made  available  to  provide  a 
more  comprehensive  evaluation  of  the  Widnall  program  by  relating  its  pre¬ 
dictions  to  available  hydrofoil  and  airfoil  data. 

Both  a  description  of  the  numerical  characteristics,  and  comparisons 
between  predicted  loading  and  experimental  results  for  the  Widnall  three- 
dimensional  calculation,  are  presented  in  this  report.  Comparisons  are 
made  for  both  steady  and  unsteady  loading  on  planar,  fully-wetted  foils. 

A  description  of  the  program  as  used,  and  instructions  for  program  use,  are 
given  in  the  appendices. 

The  results  contained  in  this  report  were  obtained  from  a  computer 
program  which  was  a  modification  of  the  original  version  received  by  the 
Center.  Changes  in  coding  were  made  by  the  authors  to  modify  input  and 
output,  to  increase  the  number  of  chordwise  integration  modes,  and  to  adapt 
the  program  to  various  computers.  Results  generated  by  the  present  program 
agree  with  previously  published  results'*'  for  equivalent  inputs.  It  is, 
therefore,  believed  that  the  present  results  are  the  same  as  would  have 
been  obtained  from  the  original  version. 

SUMMARY  OF  THE  LIFTING-SURFACE  THEORY 

The  computer  programs  described  in  this  report  are  used  to  compute 

hydrodynamic  loading  on  three-dimensional  hydrofoils  in  inviscid,  sub- 
2 

cavitating  flow  according  to  the  lifting-surface  theory  presented  by  Ashley 
et  al.l  Calculations  for  hydrofoils  in  supercavitating  flow  are  not  in¬ 
cluded.  A  brief  summary  of  the  theory  will  be  given  in  order  to  facilitate 
description  of  program  characteristics. 

The  lifting  surface  under  consideration  consists  of  a  foil  or  com¬ 
bination  of  foils  which  may  be  nonplanar.  The  surface  is  assumed  to  be  of 
zero  thickness  and  to  make  small  angles  with  the  direction  of  flow.  An 
example  of  a  single  nonplanar  foil  is  shown  in  Figure  1  along  with  a  possi¬ 
ble  coordinate  system.  The  choice  of  a  coordinate  system  is  limited  only 
by  the  requirements  that  flow  be  in  the  positive  x-direction  and  that  a 
free  surface,  when  present,  coincide  with  the  plane  z  =  0. 


The  foil  is  represented  theoretically  by  a  distribution  of  pressure 
doublets  over  the  surface  of  the  foil.  The  downwash,  or  velocity  pertur¬ 
bation  normal  to  the  undisturbed  flow,  is  a  known  quantity  determined  by 
flow  tangency  on  the  foil  and  is  given  by 

vn(x,y,z,t)  =  JJ  Ap(C,r|,C)  K(x-?,y-n, z-£, t)  dA  (1) 

A 

The  kernel  function  K  is  a  geometrically  determined  pressure  influence 
function  relating  the  pressure  jump  at  (C»H>0  to  the  velocity  induced  at 
(x,y,z).  Time  dependence  is  included  as  simple  harmonic  variation  in  the 
magnitudes  of  K  and  v  ,  using  complex  number  notation.  Therefore,  un¬ 
steady  loading  may  be  calculated  only  for  sinusoidally  varying  flow  con¬ 
ditions  or  foil  motions.  Superposition  of  Fourier  components  may  be  used 
to  treat  nonsinusoidal  time  dependence.  The  kernel  function  contains  a 
term  for  the  foil  being  treated  and  a  term  for  an  image  foil  equidistant 
from  the  free  surface  when  a  free  surface  is  present.  This  representation 
of  a  free  surface  corresponds  to  assuming  that  gravity  waves  created  by 
the  foil  are  unimportant,  i.e. ,  infinite  Froude  number. 

The  pressure  jump  distribution,  Ap,  is  found  by  inverting  Equation  (1) 
to  give  Ap  as  a  function  of  known  quantities.  First,  Ap  is  written  in  the 
form  of  a  series  of  chordwise  and  spanwise  pressure  modes  with  unknown  co¬ 
efficients  as  follows: 


Ap(9,sf)  =  8iTpU 


2  b 
c(b') 


m 


z 

L  m=0 


a  f  (b' )  cot  -z 
om  m  2 


n  m' 


+  y  y  — %-  a  f  (b')  sin  n9 
„2n  run  m 


n*l  m=0 


(2) 


The  chordwise  modes  (the  trigonometric  functions  of  0)  used  in  this  ex¬ 
pansion  were  taken  from  two-dimensional  airfoil  theory,  while  the  spanwise 


modes  f  are  left  to  be  specified  by  the  program  user  according  to  boundary 
m 

conditions  as  described  in  Appendix  B. 

The  integration  indicated  in  Equation  (1)  is  performed  numerically  in 

terms  of  the  unknown  pressure  mode  coefficients  a  .  A  set  of  linear 

nm 

equations  is  obtained  relating  the  pressure  mode  coefficients  to  the  down- 
wash  at  a  number  of  locations  on  the  foil.  The  set  of  equations  is  inverted 
to  give  the  pressure  mode  coefficients  in  terms  of  the  downwash  values. 
Finally,  the  pressure  jump  distribution  and,  therefore,  any  foil  loading 
coefficient  may  be  determined  by  inserting  the  foil's  downwash  values  into 
the  inverted  equations. 

Downwash  values  are  chosen  to  correspond  to  flow  conditions  on  the 
foil.  In  order  to  achieve  flow  tangency,  the  steady  downwash  is  chosen 
to  be  proportional  to  the  local  angle  of  attack.  The  unsteady  downwash 
consists  of  flow  components  needed  to  maintain  the  flow  in  contact  with, 
and  tangent  to,  the  foil  surface.  Calculation  of  the  downwash  is  de¬ 
scribed  in  Appendix  B. 

For  numerical  reasons,  the  downwash  must  be  specified  at  a  minimum 
number  of  points  on  the  foil  corresponding  to  the  product  of  the  numbers 
of  spanwise  and  chordwise  modes.  However,  the  downwash  may  be  specified 
at  additional  points  to  obtain  presumably  more  accurate  results  by  means 
of  a  least  square  fit  inversion  technique.  Over-specification  of  downwash 
values  is,  in  fact,  essential  to  obtaining  numerical  stability,  as  de¬ 
scribed  in  the  following  section  dealing  with  numerical  behavior. 

Final  specification  of  the  pressure  jump  distribution  on  the  foil  in 

terms  of  the  pressure  mode  coefficients  a  is  obtained  by  combining  the 

nm 

downwash  input  with  the  coefficient  matrix  obtained  by  integration  and  in¬ 
version  of  Equation  (1).  Real  and  imaginary  parts  of  the  calculation  are 
treated  separately  in  order  to  obtain  the  phase  relationship  between  pres¬ 
sure  distribution  and  displacement. 

Load  coefficients  over  various  regions  of  the  foil  are  calculated  by 
expressing  the  desired  coefficient  in  terms  of  the  pressure  jump  Ap  and 
integrating  over  the  region.  Real  and  imaginary  parts  of  the  integration 
are  treated  identically.  The  present  form  of  the  computer  program  calcu¬ 
lates  spanwise  lift,  pitch  moment,  and  roll  moment  coefficients,  chordwise 
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pressure  coefficient  distribution  for  steady  loading,  total  lift,  pitch 
moment  and  roll  moment  coefficients,  and  centers  of  pressure.  Unsteady 
load  coefficients  are  given  in  the  form  of  magnitude  and  phase  angle  by 
which  the  load  coefficient  leads  the  unsteady  displacement  vector.  Other 
load  coefficients  may  be  programmed  as  desired. 

The  numerical  solution  procedure  just  described  is  controlled  by  the 
program  user,  who  specifies  the  number  of  collocation  stations  at  which 
downwash  is  input,  the  number  of  integration  intervals,  and  the  number  of 
modes  used  in  the  pressure  representation.  Instructions  for  generating 
the  program  input  data  are  given  in  Appendix  B. 

NUMERICAL  STABILITY  OF  THE  CALCULATION 

In  order  to  discuss  the  numerical  behavior  of  the  program,  the  com¬ 
puter  code  names  of  the  numerical  variables  with  their  numerical  signifi¬ 
cance  will  be  introduced  at  this  time.  A  more  complete  description  of  the 
variables  is  given  in  Table  B.5. 

NOLT  —  Number  of  chordwise  pressure  modes  (n  in  Equation  (2)) 

NOST  —  Number  of  spanwise  pressure  modes  (m  in  Equation  (2)) 

NOCP  —  Number  of  collocation  stations 

NCP,  MP  —  Numbers  of  chordwise  integration  steps 

NI  —  Number  of  spanwise  integration  steps 

The  Widnall  program  was  exercised  with  a  wide  range  of  values  for  the 
above  quantities.  The  object  was  to  find  a  combination  of  parameter  values 
which  produced  a  numerically  stable  calculation.  The  term  numerical  sta¬ 
bility  is  used  to  indicate  a  condition  in  which  several  significant  figures 
in  the  calculated  quantities  do  not  change  when  the  input  quantities  are 
varied  by  a  small  amount. 

Four  hydrofoil  configurations  were  used  as  test  cases  for  the  numeri¬ 
cal  studies.  Rectangular  hydrofoils  of  AR  =  1,  6,  and  10  provided  most  of 
the  stability  information.  An  AR  =  6.1  foil  with  a  sweep  angle  A  of  15 
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deg  and  a  taper  ratio  T  of  0.25  represented  the  effects  of  sweep  and  taper. 
The  foils  were  given  a  NACA  a  =  1.0  mean  line  to  study  camber  loading  pre¬ 
dictions.  This  choice  of  camber  was  based  on  the  presence  of  the  NACA 
a  =  1.0  camber  line  on  most  of  the  hydrofoil  models  for  which  data  are 
available. 

To  summarize  the  results,  the  degree  of  numerical  stability  that  could 
be  achieved  was  quite  limited,  being  more  than  one  but  less  than  three 
significant  figures.  More  precisely,  the  calculations  for  the  AR  =  10 
hydrofoil  had  a  residual  instability  of  0.4  to  9  percent,  depending  on  the 
load  coefficient  considered.  Most  coefficients  became  more  stable  at 
lower  aspect  ratio;  at  AR  =  1  the  instability  was  about  1  percent  for  all 
coefficients.  Results  for  the  swept,  tapered  planform  were  similar,  except 
for  the  pitching  moment  slope  which  had  an  instability  of  27  percent. 
Specific  effects  of  the  numerical  parameters  and  model  configurations  will 
be  discussed  in  the  following  sections. 

NUMERICAL  STABILITY  OF  THE  STEADY  LOADING  CALCULATION 

A  large  number  of  calculations  were  performed  for  hydrofoils  under 
steady  flow  conditions.  Camber- induced  loading  was  investigated  separately 
from  loading  due  to  angle  of  attack.  Camber  loading  is  expressed  in  terms 
of  the  coefficients  C^(a  =  0)  and  C^(a  =  0).  Angle  of  attack  loading  is 
presented  as  the  coefficient  slopes  C^Ca  =  0)  and  C^a(ct  =  0)  • 

Number  of  Chordwise  Pressure  Modes,  NOLT 

The  stability  of  the  steady  loading  calculation  was  extremely  sensitive 
to  the  number  of  chordwise  pressure  modes  used.  The  calculated  coeffi¬ 
cients  are  shown  in  Figures  2  through  13. 

As  the  number  of  pressure  modes  was  initially  increased,  most  values 
of  lift  and  pitching  moment  began  to  approach  a  stable  condition  in  an 
almost  monotonic  manner,  as  shown  in  Figures  2  through  7  and  11  through  13. 

A  somewhat  different  pattern  was  followed  by  the  lift  slope,  which  oscil¬ 
lated  with  decreasing  amplitude  as  it  approached  an  apparently  stable  value; 
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Figure  2  -  The  Steady  Lift  Coefficient  CL(a  =  0)  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 
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Figure  4  -  The  Effect  of  the  Number  of  Integration  Steps  NCP,  MP,  and 
NI  and  Control  Points  NOCP  on  the  Steady  Lift  Coefficient  CL(a  =  0) 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 


Figure  5  -  The  Steady  Moment  Coefficient  CM(a  =  0)  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 
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Figure  6  -  The  Effect  of  the  Number  of  Chordwise  Integration  Steps  NCP 

and  MP  on  the  Steady  Moment  Coefficient  Cw(a  =  0)  for  a  Planar 

M 

Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  7  -  The  Effect  of  the  Number  of  Integration  Steps  NCP,  MP,  and 
NI  and  Control  Points  NOCP  on  the  Steady  Moment  Coefficient  CM(a  =  0) 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  8  -  The  Steady  Lift  Coefficient  Slope  CLa  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 
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Figure  9  -  The  Effect  of  the  Number  of  Chordwise  Integration  Steps  NCP 
and  MP  on  the  Steady  Lift  Coefficient  Slope  CLa  for  a  Planar 

Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  10  -  The  Effect  of  the  Number  of  Integration  Steps  NCP,  MP,  and 
NI  and  Control  Points  NOCP  on  the  Steady  Lift  Coefficient  Slope  C^ 
for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  11  -  The  Steady  Moment  Coefficient  Slope  CMa  as  a  Function  of 

the  Number  of  Chordwise  Modes  NOLT  for  a  Planar  Rectangular 
AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  12  -  The  Effect  of  the  Number  of  Chordwise  Integration  Steps  NCP 
and  MP  on  the  Steady  Moment  Coefficient  Slope  CMa  for  a  Planar 

Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  13  -  The  Effect  of  the  Number  of  Integration  Steps  NCP,  MP,  and 
NI  and  Control  Points  NOCP  on  the  Steady  Moment  Coefficient  Slope 
Cw  for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 


see  Figures  8  through  10.  These  differences  are  attributed  to  the  different 
combination  of  pressure  modes  required  to  represent  the  chordwise  pressure 
distribution  arising  from  the  two  types  of  downwash. 

Before  the  calculation  had  become  completely  stable,  however,  the  co¬ 
efficients  abruptly  shifted  away  from  the  approaching  stable  values  and, 
thereafter,  exhibited  gross  instability.  This  latter  instability  was  dis¬ 
covered  upon  examination  of  the  calculated  chordwise  pressure  distribution 
to  be  an  occurrence  of  the  Gibbs  phenomenon.  The  Gibbs  phenomenon  refers 
to  the  gross  inaccuracy  that  occurs  near  a  discontinuity  or  sharp  gradient 
when  a  function  is  represented  by  a  Fourier  series.  The  uniform  chordwise 
pressure  jump  distribution  of  the  NACA  a  =  1.0  mean  line  is  discontinuous 
at  both  leading  and  trailing  edges,  and  the  pressure  distribution  due  to 
angle  of  attack  approaches  infinity  at  the  leading  edge.  Chordwise  pres¬ 
sure  jump  distributions  for  several  numbers  of  pressure  modes  are  illus¬ 
trated  in  Figure  14.  The  growth  of  unrealistic  spikes  near  the  leading 
and  trailing  edges  as  the  number  of  modes  increases  is  evident. 

It  was  concluded  that  ten  chordwise  modes  (NOLT  =  10)  should  be  used. 
Since  the  calculation  had  not  yet  become  completely  stable  at  that  value, 
other  numerical  parameters  were  examined  in  an  attempt  to  attain  improved 
stability. 

Number  of  Chordwise  Integration  Stations,  NCP  and  MP 

Despite  the  origin  of  the  numerical  instability  in  the  chordwise 
pressure  distribution,  additional  integration  stations  in  the  chordwise 
direction  did  not  make  the  calculation  significantly  more  stable.  The 
principal  advantage  of  raising  NCP  and  MP  from  10  to  32  was  to  delay  the 
complete  breakdown  in  the  convergence  process;  this  effect  can  be  seen 
most  clearly  in  Figures  4,  7,  10,  and  13.  The  calculation  continued  to  be 
unstable  at  NCP  =  MP  =  32,  however,  as  is  apparent  in  the  figures  and 
Table  1.  A  further  increase  in  NCP  and  MP  to  64  appeared  to  have  a  sta¬ 
bilizing  effect  on  the  camber-induced  loading  (Figures  4  and  7)  but  a  de¬ 
stabilizing  effect  on  the  loading  slopes  (Figures  10  and  13).  These 
unsuccessful  results  with  larger  numbers  of  chordwise  integration  stations 
led  to  the  selection  of  NCP  =  MP  =  10  as  baseline  values. 
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LOCAL  LIFT  COEFFICIENT  DISTRIBUTION  A  p/q 


Figure  14  -  Numerical  Instability  of  the  Chordwise  Distribution  of  Lift 
at  the  Midspan  of  an  AR  =  10  Rectangular  Hydrofoil  at  Infinite  Depth 


Figure  14a  -  Lift  due  to  a  NACA  a  =  1.0  (C^  =  0.3)  Camber  Line 


Figure  14  (Continued) 


NOCP  =60  (20  CHORDWISE,  3  SPANWISE) 
NOST  =  2 
NCP  =MP  =  10 
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Figure  14b  -  Lift  on  an  Uncambered  Foil  at  an  Angle  of 
Attack  of  1.0  Radian 


TABLE  1  -  NUMERICAL  STABILITY  OF  CALCULATED  LOAD  COEFFICIENTS  FOR  A 
PLANAR  RECTANGULAR  HYDROFOIL  AT  INFINITE  DEPTH 


AR 

NCP.MP 

NOLT 

CL(ot-0) 

cM(a»0) 

CLa{rad_1) 

CMa(rad_1) 

Percent 

Change 

Acl 

acm 

AC. 

La 

ACMa 

10 

10 

a 

0.2423 

-0.01312 

5.331 

1.471 

- 

- 

- 

- 

D 

0.2566 

-0.01552 

4.929 

1.534 

+5.9 

-18.3 

-7.5 

+4.3 

m 

0.2526 

-0.01517 

5.497 

1.490 

-1.6 

+2.3 

+11.5 

-2.9 

8 

0.2700 

-0.01654 

5.013 

1.526 

+6.9 

-9.0 

-8.8 

+2.4 

9 

0.2656 

-0.01660 

5.541 

1.539 

-1.6 

-0.4 

+10.5 

+0.9 

10 

0.2693 

-0.01835 

5.394 

1.607 

+1.4 

-10.5 

-2.7 

+4.4 

11 

0.2681 

-0.01679 

5.498 

1.465 

-0.4 

+8.5 

+1.9 

-8.9 

12 

0.3670 

-  0.00692 

3.876 

1.297 

+36.9 

+58.8 

-29.5 

-11.5 

32 

10 

0.2797 

-0.01746 

5.082 

1.543 

- 

- 

- 

- 

11 

0.2764 

-0.01721 

5.487 

1.516 

-1.2 

+1.4 

+8.0 

-1.3 

12 

0.2895 

-0.01805 

5.112 

1.537 

+4.7 

-4.9 

-6.8 

+1.4 

13 

0.2884 

-0.01795 

5.368 

1.525 

-0.4 

+0.6 

+5.0 

-0.8 

14 

0.2995 

-0.01861 

5.087 

1.539 

+3.8 

-3.7 

+5.2 

+0.9 

6 

10 

8 

0.2139 

-0.01904 

4.234 

1.250 

- 

- 

- 

- 

9 

0.2094 

-0.01885 

4.581 

1.238 

-2.1 

+1.0 

+8.2 

-1.0 

0.2184 

-0.02024 

4.338 

1.274 

+4.3 

-7.4 

-5.3 

+2.9 

0.2150 

-0.01947 

4.566 

1.223 

-1.6 

+3.8 

+5.3 

-4.0 

B 

10 

8 

0.08276 

-0.02435 

1.454 

0.4944 

- 

- 

- 

- 

■ 

9 

0.08196 

-0.02416 

1.471 

0.4904 

-1.0 

+0.8 

+1.2 

-0.8 

0.08420 

-0.02486 

1.456 

0.4953 

+2.7 

-2.9 

-1.0 

+1.0 

■ 

0.08354 

-0.02462 

1.470 

0.4901 

-0.8 

+1.0 

+1.0 

-1.0 

NOTE: 

The  above  results  were  obtained  with  60  control  points  (NOCP-60)  arranged  as  in 
Table  2.  Also,  NOST-2  and  NI-10. 


Number  of  Collocation  Points,  NOCP 

In  general,  a  large  number  of  collocation  or  control  points  are  re¬ 
quired  in  the  chordwise  direction,  while  relatively  few  suffice  in  the  span 
wise  direction.  Exploratory  calculations  were  made  using  a  basic  array  of 
20  chordwise  and  3  spanwise  rows  of  control  points  (NOCP  =  60)  located  as 
given  in  Table  2.  Evenly  distributed  chordwise  rows  of  points  were  removed 
or  spanwise  rows  added  as  required. 
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TABLE  2  -  INPUT  PARAMETER  BASELINE  VALUES  REQUIRED 
FOR  APPROXIMATE  NUMERICAL  STABILITY  OF  THE 
STEADY  LOADING  CALCULATION 


Program 

Variable 

Input  Value 

■’*'  i 

60 

10 

NOST 

2 

NCP 

10 

HP 

10 

NI 

10 

XOC 

0.025,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4, 

0.45,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9, 

0.95,0.975 

SOS 

0.35,0.65,0.9 

Two  minimum  conditions  were  found  necessary  to  avoid  marked  insta¬ 
bilities  in  the  calculation:  the  number  of  rows  of  control  points  in  both 
chordwise  and  spanwise  directions  had  to  exceed  the  number  of  pressure  modes 
in  the  respective  directions.  This  result  is  indicated  in  Table  3.  The 
effects  of  additional  control  points  beyond  this  minimum  number  are  shown 


TABLE  3  -  EFFECT  OF  CONTROL  POINT  ARRANGEMENT  AND  NUMBER  OF  PRESSURE  MODES 
ON  THE  NUMERICAL  STABILITY  OF  STEADY  LIFT  COEFFICIENT  C  (a=0) 

Li 

(AR=10  rectangular  foil;  NACA  16-300  (a=1.0)  profile;  h/c  =  °°;  NCP= 
MP=32;  NI=10) 


Number  of 
Control 

Po lnts 
(NOCP) 

Number  of 
Chordwise  Rows 
of 

Control  Points 

Number  of 
Chordwise 
Pressure  Modes 
(NOLT) 

Number  of 
Spanwise  Rows 
of 

Control  Points 

Number  of 
Spanwise 
Pressure  Modes 
(NOST) 

Lift 

Coefficient 

CL 

60 

20 

8 

3 

4 

Unstable 

60 

20 

8 

3 

2 

Stable 

84 

7 

13 

12 

4 

Unstable 

84 

7 

10 

12 

4 

Unstable 

84 

7 

7 

12 

4 

Stable 

84 

7 

4 

12 

4 

Stable 

in  Figures  15  through  18.  Adding  rows  of  control  points  in  the  chordwise 
direction  improved  the  stability  of  the  calculation,  while  increasing  the 
number  of  spanwise  rows  of  control  points  beyond  the  minimum  had  no  effect. 
The  calculations  showed  that  20  chordwise  and  3  spanwise  rows  were  suffi¬ 
cient  to  eliminate  the  more  severe  instabilities  in  the  load  coefficients. 

Number  of  Spanwise  Pressure  Modes,  NOST 

Two  spanwise  pressure  modes  (NOST  =  2)  were  sufficient  for  good  nu¬ 
merical  stability,  as  shown  in  Figures  19  through  22.  This  result  occurred 
because  the  spanwise  foil  loading  closely  resembled  the  first,  elliptically- 
shaped  pressure  mode.  (The  modes  used  for  this  test  case  were  the  same  as 
those  in  Subroutine  FUNCTN  in  Programs  2  and  3,  Appendix  C.)  Note  that 
other  spanwise  loading  conditions  may  require  additional  spanwise  pressure 
modes  for  accurate  representation. 

Number  of  Spanwise  Integration  Stations,  NI 

This  variable  was  not  extensively  evaluated,  in  view  of  the  relative 
insensitivity  of  the  calculated  coefficients  to  the  number  of  spanwise 
pressure  modes.  Only  two  calculations  were  made  with  a  value  of  NI  other 
than  10.  In  these  calculations,  shown  in  Figures  4,  7,  10,  and  13,  NI  was 
increased  from  10  to  64  while  NCP  and  MP  were  simultaneously  increased  from 
32  to  6\.  Since  the  results  did  not  change  noticeably,  it  was  concluded 
that  10  spanwise  integration  stations  (NI  =  10)  are  adequate. 

Effect  of  Foil  Aspect  Ratio 

The  stability  of  the  calculation  improved  as  aspect  ratio  decreased. 

As  derived  from  Table  1,  calculations  made  at  aspect  ratios  of  10,  6,  and  1 
had  the  following  ranges  of  variation  when  the  value  of  N0LT  was  increased 
from  10  to  11. 


Aspect 

Ratio 


Range  of  Change 
in  Load  Coefficients 
(percent) 


(a  =  0) 


» 


Figure  15  -  The  Effect  of  the  Number  of  Chordwise  Control  Points  j  and 
Spanwise  Control  Points  k  on  the  Steady  Lift  Coefficient  C^(a  =  0) 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 


Figure  16  -  The  Effect  of  the  Number  of  Chordwise  Control  Points  j  and 

Spanwise  Control  Points  k  on  the  Steady  Moment  Coefficient  C..(oi  =  0) 

M 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 
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Figure  17  -  The  Effect  of  the  Number  of  Chordwise  Control  Points  j  and 
Spanwise  Control  Points  k  on  the  Steady  Lift  Coefficient  Slope  C 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 


Figure  18  -  The  Effect  of  the  Number  of  Chordwise  Control  Points  j  and 
Spanwise  Control  Points  k  on  the  Steady  Moment  Coefficient  Slope  C^a 

for  a  Planar  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth 


NOCP  =  84 
NOLT  =  6 

NCP  =  MP  =  Nl  =  10 


Figure  19  -  The  Steady  Lift  Coefficient  C  (a  =  0)  as  a  Function  of  the 

la 

Number  of  Spanwise  Modes  NOST  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 


NOCP  =  84 
NOLT  =  6 
NCP  =  MP  =  10 


Figure  20  -  The  Steady  Moment  Coefficient  C.,(a  =  0)  as  a  Function  of  the 

M 

Number  of  Spanwise  Modes  NOST  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 


Figure  21  -  The  Steady  Lift  Coefficient  Slope  CLa  as  a  Function  of  the 

Number  of  Spanwise  Modes  NOST  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 


Figure  22  -  The  Steady  Moment  Coefficient  Slope  CMa  as  a  Function  of  the 

Number  of  Spanwise  Modes  NOST  for  a  Planar  Rectangular  AR  =  10 
Hydrofoil  at  Infinite  Depth 
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Thus,  a  large  reduction  in  the  size  of  the  maximum  numerically  induced 
shift  occurred  as  aspect  ratio  was  reduced,  although  the  minimum  changes 
remained  at  about  1  percent.  This  improvement  can  also  be  seen  in  the 
chordwise  distribution  of  lift,  shown  in  Figure  23.  Substantially  smaller 
spikes  occurred  near  the  leading  and  trailing  edges  as  aspect  ratio  was 
reduced.  The  improved  stability  may  have  been  due  to  the  change  in  the 
shape  of  the  loading  distribution  that  results  from  the  increasingly  three- 
dimensional  flow  and  is  more  easily  represented  by  the  chordwise  pressure 
modes  used. 

The  breakdown  of  the  convergence  process  seen  at  AR  =  10  was  also 
delayed  to  higher  values  of  NOLT  at  lower  aspect  ratio.  This  result  is 
displayed  in  Figures  24  through  31.  Even  at  the  higher  values  of  NOLT, 
however,  the  calculated  coefficients  continued  to  vary  substantially  rather 
than  become  stable.  Since  the  sensitivity  to  numerical  parameters  de¬ 
creased  with  decreasing  aspect  ratio,  it  was  concluded  that  the  previously 
selected  parameter  values  were  adequate  for  aspect  ratios  of  less  than  10. 

Effect  of  Finite  Depth 

The  influence  of  shallow  submergence  on  the  numerical  representation 
was  estimated  by  a  series  of  calculations  made  for  the  aspect  ratio  10  foil 
at  a  relatively  shallow  depth  of  0.2  chord.  The  variations  of  numerical 
parameters  and  the  resulting  changes  in  loading  are  shown  in  Table  4.  The 
variations  were  similar  to  those  occurring  under  infinite  depth  conditions. 
It  is  concluded  that  numerical  stability  is  not  affected  by  depth  of 
submergence. 

Effects  of  Sweep  and  Taper 

Calculations  were  made  for  an  AR  =  6.1  foil  with  T  =  0.25,  A  =  15 
deg,  and  a  NACA  16-300,  a  =  1.0  section  in  the  streamwise  direction.  Table 
5  shows  that  the  calculated  coefficients,  with  one  exception,  changed  by 
less  than  7  percent  from  values  calculated  with  previously  recommended 
parameter  values,  as  the  numerical  parameters  varied.  However,  the  ex¬ 
ception,  CMa>  oscillated  by  up  to  27  percent  as  NOLT  varied.  This  large 
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Figure  23  -  The  Effect  of  Aspect  Ratio  on  the  Chordwise  Distribution  of 
Lift  at  the  Midspan  of  a  Rectangular  Hydrofoil  at  Infinite  Depth 


Figure  23  (Continued) 


Figure  23b  -  Lift  on  an  Uncambered  Foil  at  an  Angle  of 
Attack  of  1.0  Radian 
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Figure  24  -  The  Steady  Lift  Coefficient  C^(a  =  0)  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  6  Hydro¬ 
foil  at  Infinite  Depth 
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Figure  25  -  The  Steady  Moment  Coefficient  C^(a  -  0)  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  6  Hydro¬ 
foil  at  Infinite  Depth 
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Nl  =  10 


Figure  26  -  The  Steady  Lift  Coefficient  Slope  C  as  a  Function  of  the 

LiOt 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  6  Hydro¬ 
foil  at  Infinite  Depth 
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Figure  27  -  The  Steady  Moment  Coefficient  Slope  CMa  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  6 
Hydrofoil  at  Infinite  Depth 
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Figure  30  -  The  Steady  Lift  Coefficient  Slope  CLa  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  1 
Hydrofoil  at  Infinite  Depth 


Figure  31  -  The  Steady  Moment  Coefficient  Slope  CMa  as  a  Function  of  the 

Number  of  Chordwise  Modes  NOLT  for  a  Rectangular  AR  =  1 
Hydrofoil  at  Infinite  Depth 
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TABLE  4  -  EFFECT  OF  FOIL  SUBMERGENCE  ON  THE  NUMERICAL 
STABILITY  OF  THE  STEADY  LOADING  CALCULATION 

(AR=10  rectangular  foil;  NACA  16-300  (a=1.0)  profile;  h/c  =0.2 


Parameter 

Baseline 

Value 

New 

Value 

Change  in 

Load  Coefficient  (Percent) 

CL(a=0) 

CLa 

CM(a=0) 

CMa 

NOCP 

60 

80 

+3.3 

+2.5 

NOLT 

10 

11 

-1.0 

-6.7 

NOST 

2* 

3* 

-1.3 

E9 

0.0 

NCP=MP 

10 

32 

+3.2 

-4.8 

-3.3 

10 

64 

+3.2 

-4.8 

-3.3 

NI 

10 

32 

0.0 

0.0 

0.0 

0.0 

*NOCP= 

80,  with  four  spanwise  rows  of 

control 

points. 

TABLE  5  -  EFFECT  OF  SWEEP  AND  TAPER  ON  THE  NUMERICAL 
STABILITY  OF  THE  STEADY  LOADING  CALCULATION 

(AR=6.1;  T  =  0.25;  A  =  15  deg;  NACA  16-300  (a=1.0)  streamwise 
profile;  h/c'  =  °°) 


Parameter 

Baseline 

Value 

New 

Value 

Change  in 

Load  Coefficient  (Percent) 

CL(a=0) 

CLa 

CM(a=0) 

CMa 

NOCP 

60 

80 

-5.8 

+6.7 

NOLT 

10 

11 

B9 

+2.9 

-2  7 

NOST 

2* 

3* 

0.0 

+1.0 

+5.7 

NCP'MP 

10 

32 

0.0 

KB 

+1.0 

+4.7 

*NOCP= 

80,  with  four  spanwise  rows  of 

control 

points. 

variation  indicated  that  the  C  calculation  was  not  stable,  and  that  con- 

Mot 

siderable  computer  run  time  would  be  required  to  achieve  stability.  It  is 

concluded  that  stability  of  the  load  coefficients  CL(a  =  0),  CLa,  and 

CM(a  =  0)  is  not  affected  by  sweep  and  taper,  but  that  the  stability  of 

C  »  is  greatly  reduced. 

MOt 


NUMERICAL  STABILITY  OF  THE  UNSTEADY 
LOADING  CALCULATION 

In  an  abbreviated  numerical  survey,  it  was  found  that  the  unsteady 
calculation  was  more  sensitive  than  the  steady  calculation  to  the  number  of 
spanwise  pressure  modes  used.  In  other  respects,  the  calculation  was 
similar  to  the  steady  loading  calculation. 

The  sensitivity  to  spanwise  pressure  modes  may  have  resulted  from  the 
choice  of  the  model.  An  AR  =  10  rectangular  foil  was  used,  with  unsteady 
loading  produced  by  foil  oscillation  simulating  bending  and  torsional 
motion.  The  oscillation  amplitudes  were  symmetrical  about  the  midspan  and 
varied  linearly  from  0  at  the  midspan  to  a  maximum  at  the  foil  tips.  This 
type  of  motion  produced  greater  spanwise  variation  in  loading  than  the 
constant  camber  and  angle  of  attack  used  for  the  steady  loading  study.  A 
relatively  high  reduced  frequency  of  1.97  was  used  to  accentuate  the  un¬ 
steady  aspects  of  the  calculation. 

Since  preliminary  calculations  showed  a  sensitivity  to  the  number  of 
spanwise  pressure  modes,  a  baseline  array  of  15  chordwise  and  4  spanwise 
rows  of  control  points,  located  as  given  in  Table  6,  was  used.  The  magni¬ 
tudes  and  phase  angles  of  CT  and  were  examined  for  stability:  camber 

L  M 

lift  produces  no  unsteady  component. 

As  shown  in  Table  7,  the  calculated  unsteady  loading  showed  good  sta¬ 
bility  when  numerical  parameters  were  varied  from  the  baseline  values.  The 
requirement  for  4  spanwise  rows  of  control  points  is  evident  in  Figure  32. 
Stability  breakdown  again  occurred  at  high  values  of  NOLT,  as  shown  in 
Figure  33.  Except  for  the  rearrangement  in  control  point  location  and  the 
addition  of  a  pressure  mode  in  the  spanwise  direction,  the  numerical  pa¬ 
rameter  values  selected  for  the  unsteady  calculation  were  the  same  as  for 
the  steady  calculation.  The  values  chosen  are  given  in  Table  6. 

SUMMARY  OF  RESIDUAL  INSTABILITY 

Both  steady  and  unsteady  loading  calculations  failed  to  achieve  a 
value  which  was  stable  to  several  significant  figures.  Estimating  residual 
instability  with  a  high  degree  of  confidence  is  not  possible,  because  most 
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TABLE  6  -  INPUT  PARAMETER  BASELINE  VALUES  REQUIRED 
FOR  APPROXIMATE  NUMERICAL  STABILITY  OF 
THE  UNSTEADY  LOADING  CALCULATION 


Program 

Variable 

Input  Value 

NOCP 

60 

NOLT 

10 

NOST 

3 

NCP 

10 

MP 

10 

NI 

10 

XOC 

0.025,0.05,0.1,0.15,0.25,0. 3,0.4,0.45,0.55, 
0.6,0.7,0.75,0.85,0.95,0.975 

SOS 

0.1,0.35,0.65,0.9 

TABLE  7  -  NUMERICAL  STABILITY  OF  THE  UNSTEADY  LOADING  CALCULATION 

(AR=10  rectangular  foil;  uncambered  profile;  h/c  =  00 ;  bending  and 
torsional  oscillation) 


Parameter 


Baseline  New  Oscillation 
Value  Value  Mode 


Bending 

Torsion 


NOST 

NCP=MP 


Change 

in  Load  Coefficient 

c 

c 

L 

4> 

M 

4> 

(percent) 

(deg) 

(percent ) 

(deg) 

4*  Bending 


Bending 

Torsion 

Bending 

Torsion 


*NOCP=75,  with  5  spanwise  rows  of  control  points. 
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Figure  32  -  The  Effect  of  the  Number  of  Spanwise  Control  Points  k  on  the 
Lift  Coefficient  Magnitude  |  CLcJ  and  Phase  <f>L  and  Moment  Coefficient 

Magnitude  |c  |  and  Phase  <(>^  for  Bending  and  Torsion  Oscillations 

for  a  Rectangular  AR  =  10  Hydrofoil  at  Infinite  Depth  with 
k  =  1.97  (Number  of  Chordwise  Control  Points  =  15; 

NOCP  =  15  *  k) 
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NOLT 

Figure  33b  -  Moment  Coefficient  Magnitude  |c  |  and  Phase  Angle  <p 


calculated  results  oscillated  around  a  changing  local  average,  rather  than 
oscillating  around  or  asymptotically  approaching  a  constant  limiting  value. 
Attaining  stability  may  require  much  smaller  numerical  intervals,  which 
were  not  attempted  for  cost  reasons,  or  may  be  limited  to  chordwise  loading 
distributions  which  are  easily  represented  by  trigonometric  modes. 

To  provide  a  rough  approximation  of  the  residual  instability,  the  in¬ 
cremental  change  in  load  coefficient,  due  to  an  increase  in  the  number  of 
chordwise  pressure  modes  (NOLT)  from  10  to  11,  was  used.  The  instabilities 
thus  determined  are  given  in  Table  8. 

COMPUTATIONAL  PRECISION 

The  Widnall  program  was  used  on  three  different  computers:  the  IBM 
7090  and  CDC  6700  computers  at  the  Center,  and  an  IBM  360  computer  at  the 
Johns  Hopkins  University  Applied  Physics  Laboratory  (APL) .  It  was  sur¬ 
prising  to  discover  that  initially  the  program  yielded  substantially  dif¬ 
ferent  results  on  the  different  computers  although  it  was  written  in  the 
same  language  (FORTRAN  IV).  These  preliminary  calculations  were  done  in 
single  precision  (approximately  eight  significant  figures)  except  for  the 
matrix  inversion  which  was  done  in  double  precision.  Line-by-line  compari¬ 
sons  of  intermediate  program  steps  revealed  the  following  reason  for  the 
discrepancies.  The  computations  of  the  matrix  elements  relating  the  pres¬ 
sure  coefficients  to  the  downwash  involve  differences  of  terms  of  similar 
magnitude.  In  some  cases,  the  differences  were  formed  from  meaningless, 
computer-generated  digits  carried  at  the  low  end  of  the  terms  being  sub¬ 
tracted.  It  was  found  that  these  imprecise  quantities  were  associated  with 
high  order  chordwise  pressure  modes,  which  are  insignificant  in  comparison 
with  modes  of  low  order.  Unfortunately,  however,  these  erroneous  numbers 
tended  to  cluster  in  rows  and  columns  of  the  calculated  matrix.  Conse¬ 
quently,  the  determinant  of  the  matrix  was  erroneous,  and  so  was  the  cal¬ 
culated  pressure  distribution. 

As  an  example  of  the  effect  of  computer  system  characteristics  in 
such  circumstances,  the  APL  IBM  360  used  two  different  system  control  tapes 
which  gave  differing  values  for  corresponding  matrix  elements  for  the 


TABLE  8  -  APPROXIMATE  NUMERICAL  INSTABILITIES  FOR  STEADY 
AND  UNSTEADY  LOADING  COEFFICIENTS 


Type 

of 

Loading 

Foil  Configuration 

Instability 

(Percent) 

CL(a=0) 

CLa 

cM(0t“0) 

CMa 

Steady 

AR=10,  rectangular 

1 

8 

9 

9 

AR»6,  rectangular 

2 

4 

5 

4 

AR=1,  rectangular 

1 

1 

1 

1 

AR=6.1,  swept  and 
tapered 

2 

7 

3 

27 

Unsteady 

AR=1Q,  rectangular, 
bending  and  torsion 

— 

3 

3 

single  precision  version  of  the  program.  With  the  use  of  double  precision 
values  (12  precise  digits)  for  all  of  the  program  variables  and  constants, 
the  two  systems  gave  identical  results.  In  fact,  the  results  also  agreed 
with  those  from  the  IBM  7090  computer  using  double  precision  (16  precise 
digits)  and  with  those  from  the  CDC  6700  (15  precise  digits). 

No  effort  was  made  to  determine  the  relationships  between  the  number 
of  chordwise  modes  and  the  required  number  of  precise  digits.  However, 
all  of  the  steady  and  unsteady  calculations  presented  in  this  report  were 
performed  using  15  precise  digits  and  the  Center's  CDC  6700. 

The  discovery  that  the  pressure  distribution  calculated  by  the  com¬ 
puter  program  might  have  been  imprecise  caused  concern  for  the  validity  of 
the  constants  used  in  the  program.  An  examination  of  the  unsteady  program, 
which  contained  more  constants  than  did  the  steady  program,  revealed  no 
more  than  a  1-percent  change  in  the  calculated  lift  and  moment  coefficients 
when  the  least  significant  figure  in  each  of  the  constants  was  randomly 
changed  by  one  unit.  Therefore,  it  was  concluded  that  the  program  constants 
were  valid  and  did  not  need  more  significant  figures. 

COMPARISON  WITH  EXPERIMENT 

To  determine  the  accuracy  of  the  Widnall  calculation  in  predicting 
hydrofoil  loading,  loading  values  were  calculated  for  a  large  number  of 
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small-scale  hydrofoil  models  for  which  experimental  data  were  available. 

Both  steady  and  unsteady  loading  conditons  were  represented  in  these  test 
cases.  In  addition,  because  only  a  small  amount  of  data  was  available  for 
unsteady  loading  on  hydrofoils,  the  unsteady  hydrofoil  data  were  supple¬ 
mented  by  a  limited  amount  of  unsteady  data  obtained  from  airfoil  models. 
With  regard  to  the  theory,  additional  perspective  was  provided  by  parallel 
calculations  made  with  other  theoretical  approaches.  In  the  case  of  a 
lifting-surface  theory,  programmed  by  Bandler^  from  an  original  development 
by  Richardson, ^  these  calculations  were  made  by  the  authors.  Other  cal¬ 
culated  theoretical  values  were  reproduced  from  previous  publications. 

Only  very  simple  types  of  hydrofoil  and  airfoil  configurations  were 
examined  to  minimize  demands  on  the  theories.  All  hydrofoils  consisted  of 
a  planar  lifting  surface  traveling  in  a  horizontal  orientation  at  varying 
depths  beneath  a  free  surface,  or  extended  perpendicularly  to  a  large  end 
plate.  The  airfoils  were  similarly  extended  from  a  large  wall.  The  foils 
were  either  actually  or  effectively  bilaterally  symmetrical.  Both  rec¬ 
tangular  and  swept,  tapered  planforms  were  included.  Most  models  had  sup¬ 
porting  struts,  and  one  had  a  fairing  (pod)  at  the  strut-foil  junction. 

In  such  cases,  the  foil  was  considered,  in  the  theoretical  representation, 
to  include  the  area  within  the  strut  and  pod  formed  by  extending  the  leading 
and  trailing  edges  of  the  foil.  No  corrections  were  made  for  effects  of 
the  strut  and  pod  on  the  predicted  foil  loading. 

STEADY  LOADING-THEORY  AND  EXPERIMENT 

Calculated  hydrodynamic  loading  due  to  steady  motion  was  compared  with 

data  from  twenty-one  different  hydrofoil  model  configurations  and  one  air- 

g_21 

foil  model.  The  models  are  described  in  Appendix  A.  Two  models  corre¬ 
sponded  to  foil  systems  on  U.S.  Navy  hydrofoil  craft;  the  others  covered 
a  wide  range  of  characteristics.  Both  measured  and  predicted  loading 
values  are  shown  in  Figures  34  through  40.  Ranges  of  experimental  un¬ 
certainty  are  given  when  available  from  the  data  source.  The  large  un¬ 
certainties  in  lift  at  zero  incidence  will  be  commented  on  shortly.  In 
this  study,  experimental  uncertainty  refers  to  the  estimated  maximum  range 
of  error  between  the  measured  value  and  the  actual  value. 
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Figure  34  -  Total  Steady  Lift  Coefficient  as  a  Function  of  Submergence 
for  Various  Hydrofoils  with  Rectangular  Planforms  (A  =  0;  T  =  L.0) 
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Figure  35  -  Total  Steady  Lift  Slope  Coefficient  as  a  Function  of  Sub¬ 
mergence  for  Various  Hydrofoils  with  Rectangular  Planforms 

(A  =  0;  T  =  1.0) 
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Figure  35a  -  AR  =  1 
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Figure  35b  -  AR  =  2 
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Figure  35  (Continued) 
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Figure  35  (Continued) 


Figure  35g  -  AR  =  8 
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Figure  36  -  Spanwise  Distribution  of  Steady  Lift  Coefficient  for  the 
NACA  0010  Hydrofoil  (A  =  0;  T  =  1.0;  AR  =  5) 
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Figure  37  -  Total  Steady  Lift  and  Lift  Slope  Coefficients  as  Functions 
of  Submergence  for  the  One-Eighth-Scale  PCH-1  Hydrofoil  Model 
(A  =  15  Degrees;  T  =  0.25;  AR  =  6.1) 
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Figure  38  -  Total  Steady  Loading  as  a  Function  of  Submergence  for  the 
One-Twelfth-Scale  AGEH-1  Hydrofoil  Model 
(A  =  35  Degrees;  T  =  0.3;  AR  =  3) 


Figure  38a  -  Lift  and  Lift  Slope  Coefficients 
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Figure  38c  -  Pitching  Moment  and  Pitching  Moment  Slope  Coefficients 
with  Foil  Pivot  at  5.9  Percent  of  Mean  Geometric  Chord 
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Figure  39  -  Comparison  of  Total  Lift  Coefficients  on  a  Hydrofoil 
and  an  Airfoil  in  Two-Dimensional  Flow 


Coefficient  on  a  Hydrofoil  in  Two-Dimensional 
Flow  at  Finite  Depth 


The  loading  predictions  correctly  accounted  for  the  qualitative  effects 
of  camber,  angle  of  attack,  depth,  and  aspect  ratio.  These  effects  are 
well  known  to  most  readers  and  will  not  be  dwelt  upon.  Quantitatively,  the 
accuracy  was  strongly  dependent  on  aspect  ratio:  it  was  generally  true 
that  the  predictions  were  of  good  accuracy  at  low  aspect  ratio  but  too  high 
at  aspect  ratios  between  3  and  10.  Comparisons  with  the  two-dimensional 
loading  data  in  Figures  39  and  40,  corresponding  to  infinite  aspect  ratio, 
could  not  be  made  directly.  However,  by  reference  to  the  prediction  at 
AR  =  10,  it  will  be  shown  that  the  calculated  lift  would  also  be  high  at 
infinite  aspect  ratio.  An  exception  to  this  trend  occurred  at  shallow  sub¬ 
mergences,  where  camber  lift  predictions  were  too  high  even  at  low  values 
of  aspect  ratio. 

Besides  the  variation  in  prediction  accuracy  at  different  aspect 
ratios,  it  is  striking  to  note  the  variation  among  the  data  themselves. 

Large  differences  in  measured  loading  occurred  for  models  which  were  either 
similar  or  only  slightly  different.  To  illustrate  this  point,  deviations 
between  predicted  and  measured  lift  at  maximum  submergence,  which  usually 
produced  the  most  accurate  predictions,  have  been  plotted  in  Figure  41. 

When  possible,  experimental  values  were  taken  from  the  edge  of  the  error 
band  closest  to  the  theoretical  prediction,  further  minimizing  the  result¬ 
ing  deviation.  This  "map"  of  prediction  accuracy  shows  that  the  experi¬ 
mental  values  differed  by  up  to  50  percent  in  camber  lift,  and  up  to  29 
percent  in  lift  slope  for  rectangular  hydrofoils  at  similar  aspect  ratios, 
and  by  even  larger  amounts  if  the  swept  and  tapered  planforms  are  included. 
These  results  suggest  that  there  are  substantial  "hidden"  differences  in 
the  model  configurations  or  errors  in  the  data. 

There  are  several  possible  explanations  for  such  wide  scatter  in  the 
data.  These  include  effects  of  struts  and  pods,  hydrodynamic  peculiarity 
of  the  NACA  a  =  1.0  camber  line,  viscous  effects  such  as  boundary  layer 
transition  and  separation,  and  unreported  experimental  uncertainty.  These 
effects  will  be  discussed  in  the  following  sections. 


DEVIATION  OF  THEORY  FROM  EXPERIMENT.  DEVIATION  OF  THEORY  FROM  EXPERIMENT, 

lXTH  XEXP)/XEXP  (PERCENT)  (XTH  ’  XEXP)/XEXP  (PERCENT) 
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Figure  41  -  Prediction  Accuracy  of  the  Widnall  Program  as  a  Function 
of  Aspect  Ratio  for  Steady  Lift  at  Deep  Submergence 


ERROR  DUE  TO  STRUT-POD  INTERFERENCE 

All  of  the  hydrofoil  models  were  supported  by  one  or  two  struts.  The 
model  representing  the  AGEH-1  main  foil  also  had  a  pod  at  the  strut-foil 
intersection.  In  all  cases,  the  calculations  considered  the  foil  alone  as 
though  it  extended  through  the  strut  or  pod.  The  effects  of  the  strut  are 
believed  to  have  been  small  relative  to  other  sources  of  error.  There  is 
no  readily  available  estimate  of  the  effect  of  a  pod  which  does  not  rotate 
with  a  foil.  It  is  apparent  that  more  consideration  may  have  to  be  given  to 
an  analysis  of  the  flow  in  this  type  of  configuration,  particularly  at 
zero  incidence. 

ERROR  DUE  TO  CAMBER  VARIATION 

Variations  in  camber  lift  performance  would  be  expected  when  NACA  a  = 

1.0  and  other  camber  lines  are  compared.  The  a  =  1.0  camber  line  has  been 
22  23 

found  ’  to  produce  only  74  percent  of  the  predicted  loading  for  an  air¬ 
foil  in  two-dimensional  flow.  This  contrasts  with  NACA  230-series  airfoil 
sections,  which  produced  108  percent  of  the  predicted  loading,  and  NACA 
four-digit  airfoil  sections,  which  produced  93  percent  of  the  predicted 
loading.  Thus,  comparisons  between  the  a  =  1.0  and  other  types  of  camber 
should  result  in  the  a  =  1.0  sections  having  about  25  percent  less  lift, 
relative  to  predicted  lift. 

Only  three  of  the  present  experimental  models  had  other  than  NACA  a  = 

11-13 

1.0  mean  lines.  Those  models  showed  about  the  same  agreement  with 

the  Widnall  predictions  of  camber  lift  as  the  other  models,  as  shown  in 

11  12 

Figures  34e  through  34g.  Two  of  the  models  ’  had  nonstandard  profiles 
about  which  no  information  was  available  on  their  expected  performance. 

The  third  model  had  an  NACA  23012  camber  line;  nonetheless,  its  camber  lift 
ranged  from  10  to  37  percent  below  the  predicted  value  (Figure  34g).  These 
results  suggest  that  there  was  no  strong  distinction  in  accuracy  among 
different  camber  lines,  but  there  were  insufficient  data  to  establish  a 
clear  trend  of  this  sort. 
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ERROR  DUE  TO  VISCOUS  EFFECTS 

The  accuracy  of  the  loading  predictions  was  not  systematically  af¬ 
fected  by  viscous  behavior  of  the  fluid  medium  (water  or  air) .  Experimental 
values  of  Reynolds  number  R^  ranged  from  2.7  *  10^  to  4.4  x  10^.  These 
values  encompassed  the  region  in  which  boundary  layer  transition  would  be 
expected.  No  correlation  was  found  between  prediction  accuracy  and  R^. 

A  recent  investigation  by  Jones*  has  uncovered  one  type  of  viscous 
effect  which  can  influence  model  loading  measurements.  Wind  tunnel  measure¬ 
ments,  shown  in  Figure  39,  displayed  a  nonlinearity  in  lift  at  low  angles 
of  attack.  Flow  visualization  observations  showed  that  the  shift  was 
caused  by  movement  of  transition  from  midchord  to  near  the  leading  edge  on 
the  suction  side  of  the  airfoil,  as  the  angle  of  attack  increased.  The 
shift  was  eliminated  when  turbulence-stimulating  strips  were  attached  near 
the  leading  edge.  Full-scale  systems  operating  at  much  higher  R^  would 
correspond  to  the  forward  transition  location. 

Jones'  observations  are  important  because  they  define  a  mechanism  by 
which  the  camber  lift  can  differ  by  as  much  as  22  percent  between  apparent¬ 
ly  similar  foil  models.  Small  imperfections  in  models,  randomly  dis¬ 
tributed,  would  trigger  transition  at  the  forward  location  in  some  instances 
and  not  in  others.  It  does  not  appear  that  the  lift  slope  would  be  af¬ 
fected  significantly,  because  the  nonlinearity  was  confined  to  a  small 
range  of  angle  of  attack. 

ERROR  DUE  TO  UNREPORTED  EXPERIMENTAL 
UNCERTAINTY 

Most  of  the  data  scatter  appears  to  be  caused  by  relatively  large  ex¬ 
perimental  uncertainty  in  the  data.  Two  model  parameters  are  especially 
suspect:  angle  of  attack  and  profile  geometry.  Angle  of  attack  is  not 
easy  to  determine,  and,  in  most  cases,  the  lift  due  to  an  angle  of  attack 
of  one  or  two  degrees  is  equal  to  the  entire  lift  due  to  camber.  Lift 
slope  values  also  would  be  affected  if  changes  in  angle  could  not  be 
measured  much  more  accurately  than  the  angle  itself. 


*Informal  communication  (Jun  1978)  with  E.A.  Jones  of  the  Canadian 
Defence  Research  Establishment  Atlantic. 
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Experimental  uncertainty  in  angle  of  attack  was  rarely  mentioned  in 
the  reported  data.  In  one  case,  a  minimum  uncertainty  of  +0.5  deg  could 

g 

be  inferred  from  the  Feldman  data.  This  uncertainty  led  to  relatively 
wide  error  bands  which  could  explain  the  decrease  in  camber  lift  as  aspect 
ratio  increased  from  6  to  8;  see  Figures  34g  and  34h.  The  tendency  for 
camber  lift  data  to  have  greater  scatter  than  lift  slope  data  also  suggests 
angle  of  attack  uncertainty,  because  slope  would  be  unaffected  by  systematic 
errors  in  angle.  These  observations  indicate  the  need  for  improvement  in 
the  quality  of  angle  of  attack  measurements.  It  would  also  be  helpful  for 
values  of  experimental  uncertainty  to  be  reported  in  all  experiments. 

Errors  in  model  profile,  originating  in  manufacturing  defects  or 
accumulations  of  dirt  during  testing,  would  affect  camber  lift  rather  than 
lift  slope.  Either  the  mean  line  itself,  or  the  transition  behavior,  could 
be  affected. 

Jones*  described  a  tendency  for  manufacturing  errors  to  remain  constant 
in  absolute  size,  and  thus  to  become  more  important  as  model  size  de¬ 
creases.  The  smaller  hydrofoil  models,  which  had  chords  as  small  as  0.10 

m  (3.9  in.),  would  seem  to  have  been  susceptible  to  profile  errors. 

19 

Offset  discrepancies  of  10  percent  were  reported  by  Coder,  although  these 
were  not  definitely  ascribed  to  manufacturing  error,  on  a  foil  with  c'  = 
0.220  m  (8.65  in.).  Again,  it  is  not  common  practice  to  report  profile 
accuracy  of  models,  and  such  information  could  be  significant. 

STATISTICAL  ANALYSIS  OF  DATA 

Despite  the  amount  of  data  scatter,  the  assembled  body  of  data  has 

24 

value  as  a  collection  of  multisample  data.  Such  data,  taken  by  indepen¬ 
dent  investigators,  are  not  likely  to  have  been  affected  by  systematic 
errors.  Given  enough  points,  the  data  will  form  a  normal  distribution  for 
which  the  mean  is  the  most  probable  value. 

Unfortunately,  there  are  not  enough  values  to  perform  such  an  analysis 
at  any  one  aspect  ratio  with  the  present  data.  Only  two  foil  configu¬ 
rations  have  had  several  independent  measurements  made  of  loading  charac¬ 
teristics.  The  AR  =  6  rectangular  hydrofoil  has  had  four  measurements 


*E.  A.  Jones,  informal  communication. 
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made  of  lift  slope;  the  measurements  for  the  AR  =  5.84  hydrofoil  could  be 

considered  a  fifth  value.  These  measurements  are  too  few  in  number  to  form 

a  normal  distribution.  It  may  be  useful,  however,  to  know  that  the  average 

of  the  five  values  for  C,  is  3.6  at  h/c  =  2.0.  The  predicted  value  is 

La 

4.15,  or  15  percent  high.  Closer  agreement  is  obtained  for  the  three  values 
corresponding  to  deeper  submergence. 

A  second  example  is  the  AGEH-1  model,  for  which  four  measurements 
of  lift  and  moment  have  been  made. ^ ^  ^  Camber  lift  and  lift  slope  pre¬ 
dictions  at  h/c  =  1.156  were  130  and  23  percent,  respectively,  above  the 
average  of  the  data  other  than  the  camber  lift  data  of  Reference  18.  The 
preponderance  of  the  data,  independent  except  that  the  same  model  was  used 
(except  for  the  addition  of  a  sealed  flap  to  the  foil  in  Reference  19) , 
strongly  implies  that  the  higher  camber  lift  value  of  Reference  1Q  was 
erroneous.  These  results,  incidentally,  reemphasize  the  less  accurate 
character  of  the  camber  lift  prediction. 

A  similar  analysis  can  be  made  of  the  two-dimensional  data,  for  which 
four  independent  values  of  lift  and  lift  slope  are  available: 


2-D 

Data 

Source 

Average 

CL(a=0) 

Average 

CLa 

25 

Lindsey  (airfoil) 

0.24 

6.0 

Jones  (hydrofoil)* 

0.21 

5.4 

Jones  (airfoil)* 

0.21 

5.4 

Robins  (hydrofoil,  h/c  =  3)^ 

0.26 

4.9 

Data  Average 

0.23 

5.4 

In  this  comparison,  the  Jones  airfoil  data  corresponded  to  fixed  transi¬ 
tion.  The  Robins  data  are  more  questionable  than  the  others  because  it 
was  not  certain  that  a  completely  two-dimensional  flow  condition  was 

26 

achieved.  It  can  be  seen  from  expressions  presented  by  Abramson  et  al. 
that  results  at  AR  =  10  are  about  80  percent  of  the  two-dimensional  values. 


*E.  A.  Jones,  informal  communication. 


This  correction  was  used  to  extrapolate  the  Widnall  predictions  from 
AR  =  10  (see  Table  1)  to  infinite  aspect  ratio.  The  predictions  of  lift 
and  lift  slope  were  then  46  and  25  percent  high,  respectively.  These  de¬ 
viations  were  included  in  Figure  41. 

To  provide  an  estimate  of  the  accuracy  of  the  Widnall  calculation 
over  the  complete  range  of  aspect  ratio,  a  third-order  polynomial  curve  fit 
has  been  made  to  the  deviations  between  theory  and  experiment  as  a  function 
of  the  reciprocal  of  the  aspect  ratio.  Only  rectangular  foils  were  in¬ 
cluded,  using  data  from  the  maximum  submergences.  The  resulting  curve, 
replotted  in  Figure  41  in  terms  of  aspect  ratio,  indicates  that  the  devia¬ 
tion  of  the  calculation  ranges  from  0  to  +47  percent  for  camber  lift  and 
from  -15  to  +24  percent  for  lift  slope,  for  aspect  ratios  of  1  and  above. 

By  applying  a  correction  equal  to  the  value  of  the  fitted  curve  to  the 
Widnall  prediction  at  any  aspect  ratio,  the  resulting  accuracy  is  expected 
to  be  within  +15  percent  of  the  value  measured  in  a  series  of  independent 
experiments.  It  appears  that  a  higher  degree  of  accuracy  could  be  ascribed 
to  the  corrected  predictions  if  more  data  were  available  to  define  the  cor¬ 
rection  in  a  statistically  more  rigorous  manner. 

In  view  of  the  large  amount  of  data  showing  the  effect  of  varying  sub¬ 
mergence,  a  further  correction  could  be  determined  which  would  extend  the 
Widnall  predictions  to  shallow  submergences  with  little  probable  loss  of 
accuracy.  This  correction  was  not  determined  in  the  present  study. 

COMPARISON  WITH  THE  BANDLER  PROGRAM 

The  Bandler  loading  calculation, ^  included  in  Figures  34,  35,  37,  and 
38,  provided  what  appears  to  be  much  more  accurate  predictions  than  the 
Widnall  calculation.  In  particular,  the  Bandler  program  predicted  sub¬ 
stantially  lower  values  of  camber  lift  than  the  Widnall  program.  However, 
in  view  of  the  probable  influence  of  viscous  effects  in  lowering  camber 
lift  for  NACA  a  =  1.0  mean  line  sections,  the  Bandler  program  (which  also 
lid  not  consider  viscosity)  is  thought  to  have  underpredicted  camber  lift, 
v  i  r  r  hi*r  deficiency  of  the  calculation  involved  numerical  instability. 

.  i  1 1  oil  at  ions  showed  that  the  numerical  stability  of  the  program 
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was  limited  to  one  or  two  significant  figures.  A  severe  limitation  in  the 
number  of  control  points,  and  hence  pressure  modes  (only  sixteen  collocation 
points  were  allowed)  prevented  proper  exploration  of  this  aspect.  It  is 
concluded  that  the  Bandler  calculation  is  not  suitable  for  loading  predic¬ 
tion  in  its  present  form,  but  should  be  reformulated  to  ensure  numerical 
stability.  This  effort  would  be  particularly  desirable  because  the  analysis 
permits  inclusion  of  foil  thickness  and,  in  an  unusual  feature  among 
existing  theories,  predicts  the  effects  of  gravity  wave  generation  (finite 
Froude  number)  on  foil  loading. 

UNSTEADY  LOADING— THEORY  AND  EXPERIMENT 

The  unsteady  loading  predictions  of  the  Widnall  program  were  compared 

with  data  from  four  hydrofoil  and  two  airfoil  model  configura- 
15  18  27—32 

tions.  ’  ’  The  models  are  described  in  Appendix  A.  Comparisons 

were  also  made  with  other  theoretical  results,  namely,  the  Woolston  et  al. 

31  32 

lifting-surface  program,  a  quasi-steady  calculation,  and  a  two- 

18 

dimensional  version  of  the  Widnall  calculation.  The  results  are  shown  in 
Figures  42  through  49.  In  each  case,  the  calculated  results  are  shown  in 
the  form  presented  by  the  experimental  investigator.  Predictions  for  each 
of  the  models  will  be  described  in  the  following  sections. 

OSCILLATING  FOILS 

Results  for  loading  produced  by  forced  oscillation  of  foils  are  shown 
in  Figures  42  through  46.  Individual  model  results  are  discussed  in  the 
following  sections. 

NACA  66-209 (S)  Hydrofoil 

The  lift  due  to  heave  oscillation  of  an  aspect  ratio  5  rectangular 
27 

hydrofoil  is  shown  in  Figure  42.  At  the  intermediate  depth  used  in  the 
calculation,  the  predicted  lift  was  always  larger  than  the  experimental 
lift,  with  the  error  increasing  at  larger  reduced  frequency.  The  maximum 
error  was  approximately  40  percent.  The  phase  angle  was  predicted  well 
except  at  shallow  depths  where  the  calculated  value  led  the  experimental 
value  by  up  to  20  deg. 
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Figure  42a  -  Magnitude  of  Lift  Coefficient 


Figure  42b  -  Phase  Angle  of  Lift  Coefficient 

Figure  42  -  Total  Unsteady  Lift  Coefficient  due  to  Heave  Oscillation 
as  a  Function  of  Reduced  Frequency  for  the  NACA  66-209 (S) 
Hydrofoil  (A  =  0;  T  =  1.0;  AR  =  5) 


*<de9>  C£(y>/<26o/c>91.7 


$  (deg)  Cg(y)/(280/c)91  7 


Figure  43  (Continued) 


Figure  43b  -  Magnitude  and  Phase  Angle  of  Lift 
Coefficient  due  to  Bending  at  k  =  0.5 


PERCENT  OF  SEMISPAN  FROM  ROOT 

Figure  43d  -  Magnitude  and  Phase  Angle  of  Pitching 
Moment  Coefficient  due  to  Bending  at  k  =  0.5 
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Figure  43  (Continued) 


Figure  43f  -  Magnitude  and  Phase  Angle  of  Lift 
Coefficient  due  to  Torsion  at  k  =  0.5 
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Figure  43  (Continued) 


PERCENT  OF  SEMISPAN  FROM  ROOT 

Figure  43h  -  Magnitude  and  Phase  Angie  oi  Pitching  Moment 
Coefficient  due  to  Torsion  at  k  =  0.5 
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Figure  44  -  Spanwise  Distributions  of  Unsteady  Loading  due  to  Bending 
and  Torsional  Oscillations  for  the  NACA  0010  Hydrofoil 
(A  =  0;  T  =  1;  AR  =  5) 


Figure  44a  -  Magnitude  and  Phase  Angle  of  Lift 
Coefficient  due  to  Bending  at  k  =  0.6 
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Figure  44  (Continued) 
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Figure  44e  -  Unsteady  Center  of  Pressure  at  k 


Figure  45  -  Total  Unsteady  Loading  due  to  Pitching  Oscillation 
as  a  Function  of  Reduced  Frequency  for  the  One-Twelfth-Scale 
AGEH-1  Hydrofoil  Model  (A  =  35  Degrees;  T  =  0.3;  AR  =  3) 
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Figure  45a  -  Magnitude  of  Lift  Coefficient,  for  Pitching 
Axis  at  32  Percent  of  the  Mean  Geometric  Chord 
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Figure  45d  -  Gain  and  Phase  Angle  of  Pitching  Moment  Coefficient, 


for  Pitching  Axis  at  32  Percent  of  the  Mean  Geometric  Chord 


Figure  45e  -  Magnitude  of  Lift  Coefficient,  for  Pitching  Axis 
at  6.2  Percent  of  the  Mean  Geometric  Chord 
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Figure  45  (Continued) 
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Figure  45f  -  Gain  and  Phase  Angle  of  Lift  Coefficient,  for 
Pitching  Axis  at  6.2  Percent  of  the  Mean  Geometric  Chord 


Figure  45g  -  Magnitude  of  Pitching  Moment  Coefficient,  for 
Pitching  Axis  at  6.2  Percent  of  the  Mean  Geometric  Chord 
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Figure  45  (Continued) 
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Figure  45h  -  Gain  and  Phase  Angle  of  Pitching  Moment  Coefficient, 
for  Pitching  Axis  at  6.2  Percent  of  the  Mean  Geometric  Chord 
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Figure  46  -  Total  Unsteady  Loading  due  to  Pitching  and  Flapping 
Oscillations  as  a  Function  of  Reduced  Frequency  for  the 
NACA  65A010  Airfoil  (A  =  0;  T  =  1;  AR  =  2) 


Figure  46  (Continued) 


Figure  46b  -  Magnitude  and  Phase  Angle  of  Pitching  Moment 
Coefficient  due  to  Pitching  about  Midchord 
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Figure  46e  -  Magnitude  and  Phase  Angle  of  Roll 
Moment  Coefficient  due  to  Flapping 


Figure  47  -  Total  Unsteady  Loading  in  a  Head  Sea  as  a  Function  of 
Wavelength  for  the  One-Twelfth-Scale  AGEH-1  Hydrofoil  Model 
(A  =  35  Degrees;  x  =  0.3;  AR  =  3) 
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Figure  47  (Continued) 


Figure  47b  -  Magnitude  and  Phase  Angle  of  Lift 
Coefficient  at  a  Speed  of  7.3  m/s 
(24.0  ft/sec) 


Figure  47d  -  Magnitude  and  Phase  Angle  of  Pitching  Moment 
Coefficient  about  an  Axis  at  6.2  Percent  of  the  Mean 
Geometric  Chord,  at  a  Speed  of  7.3  m/s  (24.0  ft/sec) 


Figure  47  (Continued) 


Figure  47g  -  Magnitude  and  Phase  Angle  of  Pitching  Moment  Coefficient 
about  an  Axis  at  6.2  Percent  of  the  Mean  Geometric  Chord,  at  a 
Speed  of  10.5  m/s  (34.6  ft/sec) 
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Figure  48  -  Magnitude  and  Phase  Angle  of  the  Total  Unsteady  Lift 
in  a  Head  Sea  as  a  Function  of  Reduced  Frequency  for  the 
NACA  66-209(S)  Hydrofoil  (A  =  0;  T  =  1.0;  AR  =  5) 


NACA  0010  Hydrofoil  and  Airfoil 

These  models  were  unusual  in  that  structural  vibration  modes  were 

simulated  rather  than  rigid-body  oscillations.  At  the  low-to-moderate 

15  28  29 

values  of  reduced  frequency  for  which  data  are  available,  ’  *  predic¬ 

tions  of  lift  were  fair  to  good  while  moment  predictions  were  fair  to  poor. 
Comparisons  are  shown  in  Figures  43  and  44.  Calculated  loading  for  the 
fully-submerged  condition  has  been  included  even  when  no  data  were 
available. 

A  change  from  0.1  to  0.5  in  reduced  frequency  of  airfoil  oscillation 
had  little  effect  on  the  accuracy  of  the  loading  calculation.  The  greatest 
discrepancy  in  calculated  phase  angle  was  a  lead  of  up  to  52  deg  in  pitch¬ 
ing  moment  at  k  =  0.1.  The  lift  magnitude  ranged  from  9  percent  below  ex¬ 
periment  to  1  percent  above  experiment,  while  pitching  moment  ranged  from 
15  percent  below  experiment  to  27  percent  above. 

In  a  remarkable  example  of  experimental  consistency,  the  fully- 
submerged  hydrofoil  data  agreed  fairly  well  with  the  airfoil  data  taken  at 
a  slightly  lower  reduced  frequency  (0.6  and  0.5,  respectively).  Much 
larger  deviations  characterized  the  comparison  with  the  surface-piercing 
hydrofoil  data,  particularly  the  moment  data.  The  lift  predictions  were 
from  27  percent  below  experiment  to  36  percent  above  experiment,  with  phase 
angles  up  to  27  deg  in  error.  The  moment  predictions  ranged  from  0  to  125 
percent  above  experiment,  with  phase  angle  as  much  as  56  deg  in  error. 

Some  of  the  discrepancy  in  the  surface-piercing  data  due  to  bending  was 
caused  by  a  ventilated  cavity  which  formed  along  the  upper  trailing  edge 
of  the  bending  model.  Predictions  of  unsteady  center  of  pressure  x  ,  a 
quantity  described  in  Appendix  B,  were  from  3  percent  to  16  percent  of  the 
chord  aft  of  experimental  values,  as  shown  in  Figure  44e. 

AGEH-1  One-Twelfth-Scale  Hydrofoil 

18 

Predicted  and  experimental  loading  on  the  AGEH-1  hydrofoil  are 
shown  in  Figure  45.  The  loads  are  given  both  in  coefficient  form  and  as 
gains.  The  gain  is  a  logarithmic  form  of  load  coefficient  as  indicated  in 


the  figures.  Because  the  calculated  gains  are  given  relative  to  calculated 
steady  loading  slopes,  which  did  not  agree  with  experiment,  the  gain  does 
not  directly  indicate  the  accuracy  of  the  unsteady  loading  predictions.  In 
the  present  case,  the  calculated  lift  gain  was  reduced  ralative  to  experi¬ 
ment  by  this  inaccuracy,  while  pitching  moment  gain  was  reduced  for  one 
pitching  axis  and  increased  for  the  other. 

Lift  was  predicted  with  fair-to-good  accuracy  over  a  wide  range  of  re¬ 
duced  frequency  k.  The  prediction  was  15  percent  above  experiment  at  low  k 
and  up  to  33  percent  below  experiment  for  k  above  unity.  The  theory  in¬ 
tersected  the  data  at  approximately  k  =  0.75.  Lift  phase  angles  were  espe¬ 
cially  accurately  predicted,  with  a  deviation  of  only  a  few  degrees  in  lag. 

Pitching  moment  coefficients  were  not  as  well  predicted.  Accuracy 
ranged  from  140  percent  above  experiment  at  low  k  to  a  large,  indeterminate 
amount  above  experiment  at  high  k.  Smaller  deviations  of  the  opposite  sign 
occurred  at  low  k  for  one  of  the  pitch  axes.  It  should  be  understood  that 
the  prediction  accuracy  will,  in  general,  vary  markedly  with  axis  location. 

Predicted  phase  angles  led  the  experimental  values  by  as  much  as  20  deg. 

18 

The  two-dimensional  version  of  the  Widnall  calculation  was  much  less 
accurate  in  predicting  the  lift  gain  and  phase  than  the  three-dimensional 
calculation.  The  two-dimensional  calculation  was  also,  in  general,  less 
accurate  in  its  predictions  of  moment  magnitude,  but  predicted  moment  phase 
angle  more  accurately.  In  view  of  its  overall  inaccuracy,  use  of  the  two- 
dimensional  calculation  would  not  be  justified  if  the  three-dimensional 
version  was  available. 

NACA  65A010  Airfoil 

30 

This  low-aspect  ratio  airfoil  was  tested  both  in  a  pitching  mode 
31 

and  in  a  flapping  mode.  Comparisons  between  theory  and  experiment  are 
shown  in  Figure  46. 

The  Widnall  prediction  tended  to  be  low  for  lift,  in  good  agreement 
with  experiment  for  pitching  moment,  and  high  for  roll  moment.  All  cal¬ 
culations  were  within  15  percent  of  the  average  experimental  values.  The 
phase  angle  predictions  and  data  varied  differently  with  reduced  frequency, 
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which  resulted  in  crossover  regions  and  in  deviations  as  high  as  17  deg. 
Unsteady  centers  of  pressure  were  calculated  within  10  percent  of  the  chord 
length.  The  Widnall  calculation  was  thus  fairly  accurate  for  these 
airfoils. 

The  Widnall  calculations  fell  within  a  few  percent  of  results  from 
the  Woolston  lifting-surface  program,  and  the  two  programs  appear  to  be 
nearly  equivalent. 

WAVE-INDUCED  LOADING 

Unsteady  lift  magnitude  and  phase  due  to  waves  was  predicted  with 
accuracy  similar  to  that  obtained  for  oscillating  foils.  Lift  magnitudes 
were  up  to  43  percent  in  error,  while  predicted  phase  angles  deviated  as 
much  as  36  deg  from  experiment.  Moment  predictions  were  less  accurate  than 
lift  predictions.  The  Widnall  lifting-surface  calculation  was  found  to  be 
more  accurate  than  either  a  two-dimensional  calculation  or  quasi-steady 
theory.  Calculations  and  data  from  three  models  are  presented  in  Figures 
47  through  49,  and  are  discussed  further  in  the  following  sections. 

AGEH-1  One-Twelfth-Scale  Hydrofoil 

Lift  and  pitching  moment  data  were  obtained  for  this  foil  in  head 

18 

seas  over  a  reduced  frequency  range  of  0.1  to  0.4.  The  Widnall  two-  and 
three-dimensional  predictions  are  compared  with  experiment  in  Figure  47. 

The  prediction  of  lift  magnitude  using  the  Widnall  three-dimensional 
calculation  was  extremely  accurate,  while  the  two-dimensional  version  was 
considerably  less  accurate.  The  three-dimensional  prediction  of  magnitude 
appeared  to  be  within  experimental  uncertainty  of  the  data  although  phase 
angles  were  consistently  20  to  30  deg  in  error. 

Moment  predictions  were  poor  for  both  calculations,  with  accuracies 
ranging  from  50  percent  above  experiment  to  60  percent  below  experiment. 
Calculated  phase  angles  were  from  70  deg  above  experiment  to  50  deg  below 
experiment . 

It  is  concluded  that  the  Widnall  three-dimensional  calculation  is 
significantly  more  accurate  than  the  related  two-dimensional  calculation 
in  predicting  lift  in  waves  at  reduced  frequencies  between  0.1  and  0.4. 
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NACA  66-209 (S)  Hydrofoil 

Calculated  and  experimental  lift  differed  by  up  to  43  percent  for  an 

27 

aspect  ratio  5  rectangular  hydrofoil  in  waves,  as  shown  in  Figure  48. 
More  accurate  predictions  occurred  at  deeper  depths.  Phase  angle  predic¬ 
tions  differed  from  experiment  by  as  much  as  16  deg,  although  scatter  in 
the  data  suggested  that  the  prediction  may  have  been  within  experimental 
uncertainty. 

NACA  16-509  Hydrofoil 

Unsteady  loading  in  both  head  seas  and  following  seas  was  measured  on 
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an  aspect  ratio  6  rectangular  hydrofoil.  The  results,  corresponding  to 
reduced  frequencies  between  approximately  0.04  and  0.36,  are  shown  in 
Figure  49.  Theoretical  results  from  the  Widnall  three-dimensional  program 
and  from  quasi-steady  theory  are  also  shown. 

The  Widnall  calculation  was  within  23  percent  of  experimental  lift 
magnitude  but  as  much  as  36  deg  different  from  measured  phase  angles.  No 
significant  difference  occurred  in  the  accuracy  of  predictions  for  head 
and  following  seas.  As  would  be  expected,  the  Widnall  results  became  in¬ 
creasingly  superior  to  quasi-steady  theory  as  reduced  frequency  increased 
(corresponding  to  a  decrease  in  wavelength  A).  It  is  concluded  that  the 
Widnall  three-dimensional  calculation  should  be  used  in  preference  to  the 
quasi-steady  calculation  at  reduced  frequenceis  greater  than  0.16. 

SUMMARY  OF  UNSTEADY  LOADING  PREDICTIONS 

Just  as  for  the  steady  loading  calculation,  the  unsteady  calculation 
showed  a  wide  range  of  accuracy  relative  to  the  available  data.  The  lift 
predictions  ranged  from  about  30  percent  below  experiment  to  40  percent 
above.  Pitching  moment  predictions  were  substantially  less  accurate,  up 
to  140  percent  in  error,  although,  generally,  they  were  within  about  60 
percent  of  experiment.  Predicted  phase  angles  for  lift  were,  as  a  rule, 
within  20  degrees  of  experiment;  pitching  moment  phase  angles  varied  much 
more  widely. 
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Possibly  due  to  the  limited  amount  of  data,  no  clear  dependence  of 
accuracy  was  found  on  the  mode  of  unsteady  loading  (heave,  pitch,  or  waves) 
or  the  foil  aspect  ratio.  On  the  other  hand,  foil  depth  and  reduced  fre¬ 
quency  did  affect  the  accuracy.  The  largest  deviations  between  theory  and 
experiment  occurred  at  shallow  depth  (Figures  42  and  44)  and  at  both  high 
and  low  reduced  frequency.  When  lift  predictions  at  deep  submergence  and 
moderate  reduced  frequency  are  considered,  for  example,  the  deviations 
reduce  to  less  than  20  percent.  Pitching  moment  deviations  remain  at  a 
substantially  higher  level  under  these  conditions.  The  effects  of  shallow 
submergence  and  high  reduced  frequency  may  be  the  reason  for  the  generally 
superior  agreement  of  the  airfoil  loading  data  with  theory:  the  airfoil 
data  did  not  include  the  extremely  shallow  submergences  nor  the  high  values 
of  reduced  frequency  that  the  hydrofoil  data  did. 

It  is  concluded  that,  as  a  rough  approximation,  unsteady  lift  predic¬ 
tions  may  be  made  with  an  accuracy  of  +20  percent  for  deeply  submerged 
hydrofoils  at  reduced  frequencies  up  to  about  1.0.  The  predicted  lift  will 
probably  be  more  than  20  percent  above  experiment  at  depths  of  one  chord 
and  less,  and  more  than  20  percent  below  experiment  at  reduced  frequencies 
above  1.0.  A  more  detailed  determination  of  the  unsteady  loading  predic¬ 
tion  accuracy  of  the  Widnall  program  cannot  be  made  in  the  absence  of  more 
systematic  data. 

Although  no  empirical  correction  to  the  unsteady  predictions  will  be 
proposed,  in  view  of  the  limited  amount  of  data,  it  appears  that  a  correc¬ 
tion  could  be  generated  to  account  for  effects  of  high  reduced  frequency. 
This  correction,  in  conjunction  with  depth  correction  developed  from  both 
steady  and  unsteady  loading  data,  would  permit  predictions  to  be  extended 
to  shallow  submergences  and  higher  reduced  frequencies  with  no  loss  of 
accuracy. 

DISCUSSION 

The  primary  result  of  this  evaluation  has  been  to  show  that  the 
Widnall  program  can  be  used  to  make  hydrofoil  lift  predictions  with  a 
moderate  degree  of  accuracy  (15  to  20  percent)  except  at  very  shallow 


submergences  and,  for  unsteady  loading,  at  very  high  reduced  frequencies. 
Pitching  moment  predictions  are  much  less  accurate,  which  is  not  surprising 
in  view  of  the  dependence  of  pitching  moment  on  both  the  magnitude  and 
chordwise  distribution  of  lift. 

The  use  of  an  empirical  correction  to  the  steady  lift  predictions  is 
of  course  an  admission  that  there  is  some  degree  of  inaccuracy  inherent  in 
the  lifting-surface  calculation.  Determination  of  the  exact  amount  of  in¬ 
accuracy  has  been  prevented  both  by  residual  numerical  instabilit  in  the 
program,  principally  at  high  aspect  ratio,  and  by  a  general  lack  of  pre¬ 
cision  and  scope  in  the  body  of  hydrofoil  data.  Of  particular  significance 
for  design  of  hydrofoil  craft  is  the  extemely  limited  amount  of  data  on 
swept,  tapered  hydrofoils.  For  the  two  such  planforms  included  in  this 
evaluation,  the  predictions  and  the  program  numerical  stability  were  both 
poor.  This  apparent  deficiency  must  be  investigated  to  determine  whether 
the  predictions  or  the  model  test  results  are  in  error. 

Some  of  the  findings  concerning  the  NACA  a  =  1.0  mean  line  suggest 
that  other  choices  of  mean  line  would  be  preferable.  Experimental  results 
show  that  the  sectional  lift  on  this  mean  line  is  only  about  three-quarters 
of  the  theoretical  value,  while  other  profiles  realize  approximately  the 
full  theoretical  value.  Furthermore,  a  modal  approach  to  describing  chord- 
wise  pressure  distributions  cannot  successfully  represent  this  type  of 
camber  lift.  As  a  result,  certain  types  of  calculation  such  as  estimating 
cavitation  inception  are  impossible.  If  a  camber  line  which  produces  less 
abrupt  increases  in  pressure  near  the  leading  and  trailing  edges  can  meet 
the  design  requirements,  more  accurate  predictions  would  probably  be  avail¬ 
able  from  the  Widnall  program  for  optimizing  other  design  parameters. 

CONCLUSIONS 

1.  The  Widnall  program  can  be  used,  in  conjunction  with  an  empirical 
correction,  to  predict  hydrofoil  lift  with  a  moderate  degree  of  accuracy 
(15  to  20  percent)  except  at  very  shallow  submergences  and  very  high  re¬ 
duced  frequencies.  Further  empirical  corrections  for  these  parameter 
ranges  can  be  derived  from  the  available  data. 
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2.  Establishment  of  a  higher  prediction  accuracy  may  be  possible  but 
is  now  prevented  by  lack  of  sufficient  high-quality  hydrofoil  laoding  data 
and  by  some  residual  numerical  instability  in  the  calculation. 

3.  The  Widnall  lifting-surface  program  is  preferable  as  a  design  tool 
to  the  Bandler  program,  the  Widnall  two-dimensional  program,  and  a  quasi¬ 
steady  theory  for  reasons  of  accuracy  or  numerical  stability. 

4.  Hydrofoil  model  data  suffer  widely  from  unreported  experimental 
error.  The  principal  source  of  error  is  the  angle  of  attack  measurement. 
Some  models  may  also  experience  a  shift  in  transition  which  results  in  an 
inaccurate  camber  lift  measurement.  Model  profiles  may  also  be  incorrect 
unless  proven  otherwise  by  measurement. 

RECOMMENDATIONS 

1.  Improved-quality  measurements  of  hydrofoil  loading  should  be  made 
to  permit  more  precise  determination  of  the  accuracy  of  the  Widnall  and 
other  calculations.  Particular  care  should  be  taken  to  obtain  accurate 
camber  lift  data,  by  assuring  accurate  measurements  of  angle  of  attack, 
correct  model  profiles,  and  using  transition  stimulation  or  flow  visualiza¬ 
tion  to  avoid  wrongly  scaled  transition  behavior.  Swept,  tapered  hydro¬ 
foils  with  and  without  pods  should  be  emphasized. 

2.  The  numerical  aspects  of  the  Widnall  program  should  be  investigated 
to  eliminate  numerical  instabilities. 

3.  The  Bandler  lifting-surface  computer  program  should  be  numerically 
extended  to  allow  determination  of  its  numerical  stability  and  accuracy. 


APPENDIX  A 

DESCRIPTION  OF  EXPERIMENTAL  MODELS 

HYDROFOIL  AND  AIRFOIL  MODELS  IN  STEADY  MOTION 

Most  of  the  steady  loading  data  used  in  the  present  study  were  obtained 
from  measurements  on  hydrofoils.  In  all,  some  twenty-one  different  model 
configurations  were  represented.  The  largest  systematic  body  of  data  per- 
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tained  to  the  DTMB  Series  HF-1  family  of  rectangular  hydrofoils.  These 
foils  had  0.229-m  (9- in.)  chords  and  NACA  16-309,  a  =  1.0  profiles.  Aspect 
ratios  of  1,  2,  3,  4,  6,  and  8  were  investigated.  The  experimental  data 
were  converted  from  the  originally  published  normal  and  longitudinal  force 
coefficients  to  lift  coefficients  for  the  present  comparison.  Experimental 
uncertainty  in  the  foil  angle  of  attack  a  was  substantial,  amounting  to  at 
least  +0.5  deg  according  to  internal  comparisons  among  the  data.  Estimates 
of  uncertainty  in  C  (a  =  0)  were  obtained  by  multiplying  nominal  values  of 
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lift  slope  C  by  the  uncertainty  in  a.  The  accuracy  of  angle  of  attack 
changes  was  not  reported,  so  that  uncertainty  in  could  not  be  de¬ 

termined.  Load  measurements  are  plotted  in  Figures  34  and  35. 

A  second  group  of  data  pertained  to  three  AR  =  6  rectangular  hydro¬ 
foils,  two  with  NACA  16-509  and  16-1009,  a  =  1.0  profiles  and  one  with  a 
NACA  23012  section.1"^  The  foils  had  0.128-m  (5.0- in.)  chords.  In  the 
original  reference,  variations  in  lift  coefficient  as  a  function  of  speed 
were  presented.  Essentially  constant  values  of  lift  coefficient  occurred 
below  the  speeds  at  which  cavitation  appeared;  these  constant  values  were 
used  in  the  present  work.  Uncertainty  in  foil  angle  of  attack  was  reported 
as  +0.2  to  -0.3  deg.  Error  bands  for  C  (a  =  0)  were  calculated  using  nomi- 

nal  values  of  C,  ,  but  could  not  be  calculated  for  CT  itself  because  the 
Lor  La 

range  of  uncertainty  for  lift  was  not  given.  Data  are  plotted  in  Figures 
34g  and  35f. 

An  additional  group  of  data  included  five  rectangular  hydrofoils  of 

varying  design.  One  was  an  uncambered  plate  with  an  elliptical  nose  and 

q 

a  bevelled  trailing  edge.  It  had  a  0.180-m  (7.07-in.)  chord  and  an  aspect 
ratio  of  1.  The  error  in  the  lift  slope  data  for  this  foil  is  indicated 
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in  Figure  35a  for  a  6.10  m/s  (20  ft/sec)  condition;  the  range  of  experi¬ 
mental  uncertainty  in  CL(a  =  0)  is  included  within  the  data  symbol  in 
Figure  34a.  The  second  and  third  foils,  of  0.204-m  (8- in.)  chord  length, 
had  NACA  64^A412,  a  =  1.0  profiles  and  aspect  ratios  of  4^  and  10; ^  see 
Figures  34d,  34i,  35d,  and  35h.  Data  given  in  Reference  10  had  been  cor¬ 
rected  for  tare  loading  and  angular  deflection  of  the  dynamometer,  but  the 

range  of  uncertainty  was  not  given.  The  fourth  and  fifth  foils  had  aspect 
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ratios  of  5.0  and  5.84;  see  Figures  34e,  34f,  and  35e.  The  AR  =  5  foil 
had  a  0.100-m  (3.9-in.)  chord  and  a  NACA  66-209(S)  profile.  The  (S)  indi¬ 
cated  a  "slight"  change  in  the  lower  offsets  near  the  trailing  edge.  Only 
one  value  of  angle  of  attack  was  used,  and,  consequently,  no  values  of  CLa 
were  obtained.  The  AR  =  5.84  foil  had  a  0.120-m  (4.7-in.)  chord  and  a 
circular  segment  profile  (flat  lower  surface,  convex  upper  surface)  of  7.5- 
percent  maximum  thickness.  No  values  of  uncertainty  were  reported  for  the 
latter  two  foils. 

Data  from  segmented  hydrofoils  with  simulated  bending  and  torsional 
deflections^  were  also  used.  The  rectangular  foils  had  0.457-m  (18-in.) 
chords  and  NACA  0010  profiles.  An  effective  aspect  ratio  of  5  and  an 
effectively  infinite  depth  were  produced  by  a  large  endplate,  which  was 
parallel  to  the  free  surface  and  perpendicular  to  the  hydrofoil  span. 
Dynamometers  at  five  spanwise  locations  enabled  spanwise  distributions  of 
loading  to  be  measured.  The  two  versions  of  the  foil  were  mechanically 
deflected  to  produce  bending  and  torsional  deflections,  respectively,  in 
both  steady  and  unsteady  modes.  Error  flags  shown  in  Figure  36  correspond 
to  the  reported  +0.2  deg  uncertainty  in  angle  of  attack. 

Published  loading  data  were  available  for  only  two  swept,  tapered 

hydrofoil  models.  One  was  a  1/12-scale  model  of  the  AGEH-1  (PLAINVIEW) 

main  T-foil.^  ^  The  load-bearing  foil  on  the  model  had  a  sweep  angle 

A  =  35.2  deg,  a  taper  ratio  t  =  0.3,  and  a  mean  chord  c  =  0.220  m  (8.65 

in.).  Data  are  plotted  in  Figure  38.  The  foil  had  a  NACA  16-(0. 353)08, 
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a  =  1.0  profile  defined  normal  to  the  32-percent  chord  line,  necessi¬ 
tating  calculation  of  streamwise  section  characteristics  for  use  in  the 
computer  calculations.  The  section  offsets  were  inadvertently  calculated 


for  a  25-percent  chord  reference  line,  but  the  difference  is  believed  to  be 
negligible.  The  relatively  large  model  pod  was  not  considered  in  the  cal¬ 
culation.  The  foil  was  rotated  in  pitch  about  several  different  axes,  as 
indicated  in  the  figure. 

The  second  swept,  tapered  foil  was  a  1/8-scale  model  of  the  forward 
T-foil  on  the  PCH-1  (HIGH  POINT)  hydrofoil  craft.  This  foil  is  designated 
as  Model  A  in  Reference  16.  It  had  A  =  15  deg,  T  =  0.35,  AR=  6.1,  and 
c  =  0.150  m  (5.9  in.).  Data  are  shown  in  Figure  37.  The  profile  was  a 
NACA  16-309,  a  =  1.0  section  defined  in  the  streamwise  direction.  No  range 
of  experimental  uncertainty  was  reported  for  either  of  the  swept,  tapered 
models . 
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Four  sets  of  measurements  ’  on  a  commonly-used  hydrofoil  section, 
the  NACA  16-309  (a  =  1.0)  profile,  in  two-dimensional  flow  were  also  in¬ 
cluded  to  represent  infinite  aspect  ratio.  Two  of  the  studies  were  done 
with  hydrofoil  models,  and  two  with  airfoil  models.  One  hydrofoil  model 

had  a  0.15-m  (6-in.)  chord  length  and  was  tested  in  the  two-dimensional 
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test  section  of  a  water  tunnel.  Airfoils  of  0.127-m  (5- in.)  and  0.762-m 

(30-in.)*  chord  were  tested  in  wind  tunnels.  Experimental  uncertainties 

were  not  given.  Data  for  one  of  the  airfoils  and  the  hydrofoil  are  com- 
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pared  in  Figure  39.  The  second  hydrofoil  had  a  0.15-m  (6- in.)  chord  and 
was  towed  beneath  the  free  surface  of  a  towing  basin;  data  are  shown  in 
Figure  40.  Struts  of  0.4-m  (15.7-in.)  chord  were  placed  at  each  end  of  the 
foil  to  approximate  two-dimensional  flow. 

HYDROFOIL  AND  AIRFORIL  MODELS  IN  UNSTEADY  MOTION 
NACA  66-209(S)  Hydrofoil 

This  model  has  been  previously  described  in  the  preceding  section  as 

27 

a  source  of  steady  loading  data.  It  was  tested  in  a  rigid  heaving  mode 

at  a  model  speed  of  6.1  m/s  (20  ft/sec)  and  a  fixed  angle  of  attack  of  1.8 

deg.  Results  are  shown  in  Figure  42.  The  lift  coefficient  used  in 

Reference  27  was  not  fully  defined  and  was  here  assumed  to  be  CTc'/(2n6  ) 

L  o 

where  6  is  the  amplitude  of  the  heave  motion. 


*Also,  informal  communication  (Jun  1978)  with  E.  A.  Jones  of  the 
Canadian  Defence  Research  Establishment  Atlantic. 
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Unsteady  loading  was  also  measured  for  this  foil  in  waves,  again  at 
a  1.8-deg  angle  of  attack;  see  Figure  48.  The  lift  data  in  waves  were  in¬ 
ferred  to  have  been  reported  in  the  form  Cl/2ttcxw,  where  a  was  the  amplitude 
of  the  effective  unsteady  angle  of  attack  due  to  the  waves.  No  experimental 
uncertainty  values  were  given. 

NACA  0010  Hydrofoil 

This  category  includes  the  two  segmented  models  also  used  for  steady 

load  measurements.  Foil  motions  were  produced  by  oscillating  the  tips  of 

the  foils  in  translation  and  rotation,  resulting  in  spanwise  mode  shapes 

approximating  the  first  bending  and  first  torsion  modes  of  a  cantilever 

beam.  The  axis  of  torsional  oscillation  was  at  the  quarter-chord  line.  A 

large  end  plate  doubled  the  aspect  ratio  when  the  foil  was  fully  submerged. 

Pitching  moments  were  reported  about  an  axis  at  31  percent  of  the  chord. 
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The  models  were  tested  at  both  full  submergence  ’  and  at  65-percent 
29 

submergence.  At  65-percent  submergence,  the  foil  should  be  referred  to 
as  a  vertical  surface-piercing  strut. 

Although  these  foils  were  tested  at  several  values  of  reduced  fre¬ 
quency,  dynamometer  inaccuracy  caused  all  data  except  that  below  k  =  0.67 
to  be  unreliable.  As  a  result,  unsteady  data  are  available  only  at  k  =  0.6; 
see  Figure  44.  The  range  of  uncertainty  in  these  remaining  data  was  not 
reported. 

NACA  0010  Airfoil 

Because  no  unsteady  data  were  available  from  the  flexible  hydrofoil 

models  mentioned  above  at  low  reduced  frequencies,  data  from  two  very 
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similar,  uncambered  NACA  0010  profile  airfoils  have  been  included  in 
this  study.  The  data  were  obtained  for  both  bending  and  torsional  oscilla¬ 
tions.  Although  similar  in  size  to  the  hydrofoil  model,  the  airfoil  model 
was  oscillated  about  the  35-percent  chord  line  during  torsional  oscilla¬ 
tions,  and  moments  were  reported  about  this  axis.  Unsteady  data  at  the 
lowest  reduced  frequency,  0.1,  and  highest  reduced  frequency,  0.5,  were 


used  for  comparison  with  theory.  These  data  are  presented  in  Figure  43. 
While  no  range  of  uncertainty  was  given,  the  data  are  believed  to  be  of 
good  accuracy. 

AGEH-1  One-Twelfth-Scale  Hydrofoil 

The  third  hydrofoil  configuration  was  a  1/12-scale  model  of  the  AGEH-1 
main  strut-pod-foil  system.  This  model  was  described  in  the  preceding 
section  on  steady  loading. 
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One  type  of  unsteady  loading  was  obtained  by  oscillating  the  foil 

in  pitch  about  either  of  two  axes  at  reduced  frequencies  k  up  to  2.  The 

axes  were  located  at  5.94  and  31.94  percent  of  c',  measured  aft  of  the 
19 

leading  edge.  These  locations  correspond  to  47  and  64  percent  of  the 
root  chord  from  the  leading  edge.  In  each  case,  the  lift  and  the  pitching 
moment  at  the  axis  were  measured.  The  experimental  procedure  did  not  com¬ 
pensate  for  inertial  effects  of  the  foil  and  linkage.  Therefore,  the 
authors  weighed  and  measured  the  model,  which  was  available  and  intact, 
and  calculated  inertial  corrections  to  the  experimental  data  published  in 
Reference  18.  These  corrections  were  not  negligible,  but  did  not  greatly 
alter  the  published  data.  The  corrected  data,  shown  in  Figure  45,  were 
used  in  this  report. 

18 

Foil  loading  was  also  measured  in  head  seas  at  a  depth  of  0.25  m 

(10  in.)  for  speeds  of  7.3  and  10.5  m/s  (24.0  and  34.6  ft/sec).  The  waves 

had  0.10-,  0.15-,  and  0.25-m  (4-,  6-,  and  10-in.)  peak-to-peak  amplitudes 

£  and  3.8-  to  12.2-m  (12.5-  to  40-ft)  wavelengths  X,  although  not  all 
w 

possible  combinations  of  amplitude  and  wavelength  were  used.  The  steady 

angle  of  attack  was  1.4  deg  at  the  lower  speed  and  -1.5  deg  at  the  higher 

speed.  Measured  loading  values  are  given  in  Figure  47.  The  downwash  used 

in  the  loading  prediction  included  both  the  horizontal  and  the  vertical 

33 

components  of  orbital  velocity.  It  was  assumed  that  the  depth  of  the 
water  was  greater  than  one-half  the  wavelength.  Since  the  towing  basin 
used  for  the  experiment  was  4.9  m  (16  ft)  in  depth,  this  assumption  was 
invalid  for  wavelengths  longer  than  9.8  m  (32  ft). 

The  experimental  uncertainty  in  the  measured  loading  was  not  given  for 
this  model. 
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NACA  16-509  Hydrofoil 

The  remaining  wave-loading  data  were  obtained  from  an  AR  =  6  rectangu- 

32 

lar  foil  with  a  NACA  16-509,  a  =  1.0  profile.  The  foil  was  experimentally 
tested  at  a  depth  of  one  chord  (76  mm  or  0.25  ft)  and  an  angle  of  attack  of 
4  deg  in  both  head  seas  and  following  seas.  Computer  predictions  for  the 
lift  were  compared  with  experimental  results  for  a  model  speed  of  3.1  m/s 
(10  ft/sec);  see  Figure  49.  The  computer  predictions  were  made  only  for 
X  =  1.5  and  2.4  m  (5  and  8  ft)  in  head  seas  and  for  X  =  0.9  m  (3  ft)  in 
following  seas.  The  values  of  k  were  calculated  using  a  calculated  wave 
propagation  speed.  In  head  seas,  k  varied  from  0.365  to  0.161  as  X  in¬ 
creased;  in  following  seas,  the  variation  was  from  0.159  to  0.035.  The 
plot  of  lift  phase  versus  wavelength  for  the  following  seas  which  appears 
in  Reference  32  is  inconsistent  in  the  labeling  of  the  ordinate  ip.  Thus 
the  experimental  phase  is  unknown  within  40  deg.  The  labeling  used  in 
this  report  assumes  that  the  phase  due  to  following  seas  is  approximately 
180  deg  greater  than  that  due  to  head  seas.  Since  the  phase  associated 
with  head  seas  was  correctly  plotted,  it  was  determined  that  the  experi¬ 
mental  phase  due  to  following  seas  clustered  around  -100  deg  and  not 
around  -60  deg.  Dynamometer  response  limitations  produced  an  uncertainty 
of  +5  percent  in  lift  amplitude  and  +5  deg  in  phase  angle. 

NACA  65A010  Rectangular  Airfoil 

The  second  airfoil  model  had  an  AR  =  2  rectangular  planform  with  a 
30  31 

NACA  65A010  profile.  ’  The  chord  length  of  the  foil  was  0.305  m  (12 

30 

in.).  A  rigid  semispan  model  of  the  foil  was  oscillated  in  pitching  and 
31 

flapping  modes  at  Mach  numbers  from  0.2  to  0.7.  The  axis  for  the  flapping 
oscillations  was  a  streamwise  line  0.152  m  (6.0  in.)  behind  the  wall  used 
for  the  image  plane;  the  pitching  axis  was  the  midchord  line.  The  mean 
angle  of  attack  was  zero.  The  experimental  and  predicted  lift  and  pitching 
moment  at  midchord  for  both  the  pitching  and  flapping  oscillations  and  the 
roll  moment  at  the  flapping  axis  for  the  flapping  oscillation  are  compared 
from  k  =  0.2  to  k  =  0.8.  The  exact  definition  of  the  experimental  phase 
for  the  pitching  oscillations  could  not  be  determined  because  of  the 
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ambiguous  notation  used  in  Reference  30,  but  was  inferred  by  comparison 
between  theory  and  experiment.  The  loading  was  measured  with  an  uncertainty 
of  +4  percent  in  magnitude  and  +3  deg  in  phase  angle. 
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APPENDIX  B 

COMPUTER  PROGRAM  USER'S  MANUAL 


The  following  material  constitutes  a  user's  instruction  manual  for 
the  Widnall  three-dimensional  lifting-surface  computer  program.  A  com¬ 
plete  description  is  given  of  the  input  and  output  procedures  and  of  all 
output  quantities.  In  subsequent  appendices,  listings  of  program  decks 
and  sample  input  and  output  are  provided. 

In  order  to  guide  a  potential  user,  it  is  noted  that  the  lengthy  in¬ 
structions  given  in  this  appendix  might  require  several  weeks  for  a  person 
who  is  already  familiar  with  FORTRAN  programming  to  learn.  The  appropriate 
sample  calculations  should  be  run  before  generating  any  original  results. 

PROGRAM  USAGE 

The  lifting-surface  theory  is  contained  in  three  computer  programs: 
the  steady  loading  program  YKWS,  the  unsteady  program  YKWU,  and  the 
auxiliary  program  YKWA.  The  programs  are  related  to  each  other  as  shown  in 
Figure  B.l.  Both  the  steady  (YKWS)  and  the  unsteady  (YKWU)  programs  cal¬ 
culate  load  coefficients  and  also  generate,  at  the  option  of  the  user, 
punched  output  containing  the  pressure  matrices  which  characterize  the 
foil  configurations  being  treated.  This  punched  output  can  then  be  used  to 
calculate  loading  for  additional  sets  of  downwash  conditions  in  relatively 
short  run  time  by  using  the  auxiliary  program  (YKWA). 

The  three  programs  are  each  composed  of  a  main  program  and  several 
subprograms.  All  programming  has  been  done  in  the  FORTRAN  Extended  Lan¬ 
guage  (Version  3.0).  The  FORTRAN  names  and  mathematical  functions  of  all 
main  programs  and  subprograms  are  shown  in  Tables  B.l  through  B.3. 

Listings  of  the  FORTRAN  coding  for  all  programs  are  given  in  Appendix  C. 

PROGRAM  INPUT  DESCRIPTION 

Input  to  the  computer  programs  consists  of  coding  in  subprograms  in 
addition  to  data  cards.  The  subprograms  specify  the  foil  geometry  and  the 
spanwise  pressure  modes.  The  data  cards  contain  foil  and  flow  parameters. 


Figure  B.l  -  Relationship  among  Programs  YKWS,  YKWU,  and  YKWA 


TABLE  B.l  -  COMPOSITION  OF  DECK  FOR  PROGRAM  1  (YKWS) 


Program 

Function 

Program  YKWS 

Main  program  (steady) 

Subroutine  ANSWER 

Calculate  steady  load  coefficients 

Subroutine  CHORD 

Calculate  chordwise  pressure  coefficient  distribution 

Subroutine  MNGLR 

Singular  spanwise  integration 

Subroutine  KERNL 

Calculate  steady  kernel  function 

Subroutine  MATINS 

Matrix  inversion 

Subroutine  CHDWS 

Chordwise  integration  (steady) 

Function  GN 

Gaussian  integration  intervals 

Function  WN 

Gaussian  weighting  factors 

Subroutine  STRUCT 

Foil  geometry 

Function  XLE 

Foil  geometry 

Function  B 

Foil  geometry 

Subroutine  FUNCTN 

Spanwise  pressure  modes 

TABLE  B. 2  -  COMPOSITION  OF  DECK  FOR  PROGRAM  2  (YKWU) 


Program 

Function 

Program  YKWU 

Main  program  (unsteady) 

Subroutine  ANSWU 

Calculate  unsteady  load  coefficients 

Subroutine  KERNLU 

Calculate  unsteady  kernel  function 

Subroutine  MNGLRU 

Singular  spanwise  integration 

Subroutine  CHDWU 

Chordwise  integration  (unsteady) 

Function  GN 

Gaussian  integration  intervals 

Function  WN 

Gaussian  weighting  factors 

Subroutine  CMPINV 

Complex  matrix  inversion 

Subroutine  MATINV 

Matrix  inversion 

Subroutine  HINT 

Auxiliary  in  unsteady  kernel  function  calculation 

Subroutine  W1NT 

Auxiliary  in  unsteady  kernel  function  calculation 

Subroutine  BIMLI 

Auxiliary  in  unsteady  kernel  function  calculation 

Function  BIML2 

Auxiliary  in  unsteady  kernel  function  calculation 

Subroutine  BESK 

Calculates  K  Bessel  function 

Subroutine  STRUCT 

Foil  geometry 

Function  XLE 

Foil  geometry 

Function  B 

Foil  geometry 

Subroutine  FUNCTN 

Spanwise  pressure  modes 
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TABLE  B. 3  -  COMPOSITION  OF  DECK  FOR 
PROGRAM  3  (YKWA) 


Program 

Function 

Program  YKWA 

Subroutine  FUNCTN 

Function  XLE 

Function  B 

Main  program  (auxiliary) 

Spanwise  pressure  modes 

Foil  geometry 

Foil  geometry 

internal  program  constants,  control  point  locations,  and  downwash  values  on 
the  foil.  The  method  of  formulating  the  required  subprograms  and  data 
values  is  described  in  the  following  sections. 

Representation  of  Lifting  Surface 

Subprogram  coding  is  used  to  supply  most  of  the  information  about  foil 
geometry  required  by  the  lifting-surface  calculation.  Any  combination  of 
foils  of  arbitrary  planform,  including  both  connected  and  unconnected 
lifting  surfaces,  can  be  treated  by  the  program.  However,  it  should  be  re¬ 
membered  that,  theoretically,  all  foil  surfaces  make  infinitesimally  small 
angles  with  the  direction  of  flow. 

The  foil  being  specified  is  located  in  the  (x,y,z)  coordinate  system 
of  Figure  1.  Flow  is  in  the  positive  x  direction.  Three  user-coded  sub¬ 
programs  are  used  to  specify  the  foil  coordinates  in  terms  of  two  new 
parameters,  the  section  number  and  the  position  along  the  span  of  the  foil, 
given  as  fraction  of  span.  Before  coding  for  these  subprograms  can  be 
written,  however,  integration  regions  on  the  foil  and  section  numbers  must 
be  chosen  to  be  consistant  with  each  other  and  also  with  the  fraction  of 
span  coordinate  directions.  Each  of  these  quantities  will  be  explained 
separately.  Names  of  variables  given  in  capital  letters  refer  to  computer 
input  quantities  which  are  listed  in  Table  B.4. 

Section  Number,  N.  The  foil  is  considered  to  be  composed  of  N  sections  of 
arbitrary  extent.  For  convenience,  N  should  be  small  and  the  sections 


TABLE  B. 4  -  FORTRAN  INPUT  NOTATION 


Input  Definition 

AF(I)  Spanvlse  length  of  1th  section,  in  units  of  reference 

semichord 

ARM  X-coordinate  about  which  pitching  moment  is  calculated 

(zero  if  left  blank) 

B(N,S)  Length  of  foil  semichord  at  the  fractional  spanwlse  position 

S  on  foil  section  N  (defined  in  FUNCTN  B(N,S)) 

EM  Mach  number;  zero  f  hydrofoils 

ETA  Half-width  of  the  singular  region  (Region  I);  minimum  of  0.1, 

SOS,  and  1.0  -  SOS 

F(I)  Value  of  the  Ith  spanwlse  pressure  mode  (defined  in  Subroutine 

FUNCTN  (N,S,F)  for  fractional  spanwlse  position  S  on  foil 
section  N) 

1PUNCH  Punched  output  indicator:  IPUNCH  »  0  generates  no  punched 

output;  IPUNCH  ■  1  generates  punched  output  required  for 
running  Program  YKWA 

J3  Free  surface  indicator:  J3  ■  0  signifies  infinite  depth; 

J3  ■  1  signifies  finite  depth  (a  value  for  the  finite  foil 
depth  must  be  given  in  Subroutine  STRUCT) 

MP  Number  of  chordwlse  integration  steps  in  the  singular  region 

(may  be  6,  8,  10,  16,  or  32) 

N  Section  number  (N  <_  6) 

NCASE  Number  of  foil  configurations 

NCMBR  Number  of  downwash  sets  for  each  foil  configuration  (given 

once  but  applies  to  all  configurations);  NCMBR  -  0  is  valid 

NCP  Number  of  chordwlse  Integration  steps  in  the  nonsingular 

regions  (may  be  6,  8,  10,  16,  or  32) 

NI(I)  Number  of  spanwlse  integration  steps  in  the  Ith  region  (may 

be  0,  6,  8,  10,  16,  or  32);  zero  prevents  integration  over 
regions  which  are  not  needed  to  represent  foil 

NNNN  Running  index  used  to  consecutively  number  the  NCASE  foil 

configurations 

NOCP  Total  number  of  control  points  over  all  sections  (NOCP  NOST 

*  NOLT;  maximum  values  are  180  for  steady  loading  and  90  for 
unsteady  loading) 

NOLT  Number  of  internally  generated  chordwlse  pressure  modes 

(NOLT  <  20) 

NOST  Number  of  spanwlse  pressure  modes  (NOST  £  4)  specified  in 

Subroutine  FUNCTN 

NS  Number  of  foil  sections  (NS  6) 

PSI  Arc  tangent  of  the  slope  dz/dy  of  the  foil  surface  in  the 

y-z  plane  when  looking  upstream  (see  Figure  B.3) 

RF  Reduced  frequency;  buj/U 

SF(I)  Ratio  of  the  span  of  the  Ith  foil  section  to  the  sum  of  the 

spans  of  all  foil  sections 

SOS  Spanwlse  location  of  control  point  on  a  given  foil  section, 

specified  as  fraction  of  span  from  the  end  of  the  section 
having  S  *  0 

SPAN  Sum  of  AF  values,  equal  to  the  total  foil  span 

VR(I)  Real  part  of  the  downwash  at  the  Ith  control  point 

VI (I)  Imaginary  part  of  the  downwash  at  the  Ith  control  point 

XOC  Chordwlse  location  of  control  point,  specified  as  fraction  of 

local  chord  aft  of  leading  edge 

YS  Fraction  of  span  at  which  chordwlse  pressure  jump  distribution 

is  to  be  calculated 
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should  correspond  to  areas  of  simple  geometry  on  the  foil  such  as  semispans, 
supporting  struts,  etc.  When  geometrical  symmetry  between  two  sections  of 
a  foil  exists,  however,  input  simplification  can  be  achieved  by  extending 
the  calculation  over  the  similar  section  without  giving  it  a  section  number. 
This  procedure  involves  manipulation  of  integration  regions. 

Fraction  of  Span,  S.  Each  numbered  section  of  the  foil  must  be  given  a 
signed  coordinate  S  which  specifies  the  fractional  spanwise  distance  b' 
from  one  end  of  the  section  to  the  other.  The  span  S  may  originate  at 
either  end  and  will  range  from  0  to  1  in  value.  The  computer  calculation 
assumes  that  negative  values  of  S  correspond  to  another,  unnumbered,  sec¬ 
tion  which  the  program  user  may  wish  to  have  included  as  part  of  the  com¬ 
plete  foil.  This  unnumbered  section  would  have  a  geometry  describable  by 
the  coding  used  for  the  numbered  section  with  negative  S  substituted  for 
positive  S. 

Integration  Regions.  The  entire  foil  configuration  is  divided  into  four  or 
more  regions  for  spanwise  integration  of  the  kernel  function.  The  regions 
change  in  location  for  each  control  point,  and,  therefore,  can  be  visual¬ 
ized  for  only  one  control  point  at  a  time.  A  narrow  chordwise  strip  cen¬ 
tered  at  the  control  point  is  called  Region  I.  ETA  is  the  spanwise 
half-width  of  Region  I.  The  remainder  of  the  foil  section  outboard  (larger 
S)  of  Region  I  is  Region  II,  while  the  inboard  part  is  Region  IV.  Region 
III  encompasses  the  unnumbered,  negative-S  counterpart  (of  the  numbered 
wing  section  composed  of  Regions  I,  II,  and  IV).  Other  even-numbered 
regions  consist  of  the  remaining  foil  sections,  while  odd-numbered  regions 
are  the  negative-S  counterparts  of  the  remaining  sections.  A  typical  set 
of  integration  regions  is  shown  on  the  foil  in  Figure  B.2. 

Integration  regions  are  designed  to  permit  special  integration  tech¬ 
niques  to  be  used  in  the  vicinity  of  each  control  point,  where  a  singu¬ 
larity  occurs  in  the  kernel  function.  Because  of  the  way  in  which  the  in¬ 
tegration  regions  are  chosen,  the  data  inputs  NI(J),  which  designate  the 
number  of  spanwise  integration  stations  in  the  Jth  region,  are  also  used  to 


indicate  whether  a  region  exists.  A  zero  value  for  all  NI(J)  where  J  is 
odd  and  greater  than  1  prevents  integration  over  negative-S  counterparts 
of  specified  foil  sections.  Because  the  integration  regions  cover  dif¬ 
ferent  sections  for  different  control  points,  all  or  none  of  the  sections 
must  have  negative-S  counterparts.  Therefore,  either  all  values  of  NI(J) 
must  be  nonzero,  or  all  N1(J)  for  J  odd  and  greater  than  1  must  be  zero. 

After  a  foil  configuration  has  been  divided  into  sections  and  inte¬ 
gration  regions  have  been  selected  to  comprise  the  full  extent  of  the  foil, 
the  required  subprograms  must  be  written  to  describe  foil  geometry  and 
spanwlse  pressure  modes.  In  describing  the  foil  geometry,  all  foil  dimen¬ 
sions  are  given  in  units  of  a  reference  length,  which  is  usually  taken  as 
the  root  semichord  of  one  of  the  foil  sections.  It,  therefore,  is  con¬ 
venient  to  place  the  foil  in  the  (x,y,z)  coordinate  system  shown  in  Figure 
1  with  the  reference  semichord  on  the  x-axis  in  the  interval  -1  x  0, 
although  such  a  choice  of  location  is  not  necessary.  Three  subprograms 
describe  the  location  and  planform  of  the  foil  in  this  frame  of  reference. 
The  subprograms  and  their  variables  are  as  follows. 

Subroutine  STRUCT.  Subroutine  STRUCT  (N, S, Z,Y,PSI)  defines  the  y  and  z 
coordinates  and  the  slope  at  each  point  on  the  foil.  The  parameter  PSI  is 
defined  to  be  the  arc  tangent  of  the  slope  dz/dy  of  the  foil  surface  in 
the  y-z  plane  when  looking  upstream  as  shown  in  Figure  B.3;  N  denotes  the 
foil  section;  Y  and  Z  represent  the  y  and  z  coordinates  of  the  foil.  The 
fractional  distance  along  the  span  of  the  section,  S,  denoted  as  b1  in 
Equation  (2),  can  be  chosen  to  originate  at  either  end  of  the  section.  The 
fractional  distance  S  will  be  set  equal  to  -S  when  unnumbered  counterparts 
to  section  N  are  being  described.  If  more  than  one  foil  configuration  is 
treated  with  a  single  computer  run,  the  variable  NNNN,  placed  in  labeled 
common  blocks  CNRA,  provides  a  running  index  to  number  the  various 
conf igurations . 

Foils  may  be  treated  either  with  or  without  a  free  surface.  The  J3 
input  described  below  indicates  whether  or  not  a  free  surface  is  to  be 
included.  The  free  surface  is  represented  by  an  image  foil  which  is 
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NOTE:  INTEGRATION  REGIONS  DENOTED 
BY  ROMAN  NUMERALS 


Figure  B.2  -  Integration  Regions  for  a  Typical  Control  Point 


Figure  B.3  -  View  of  Nonplanar  Foil  Looking  Upstream 


generated  by  using  the  given  foil  geometry  with  z-values  of  opposite  sign. 
Therefore,  foils  beneath  a  free  surface  must  be  given  z-coordinates  so  that 
the  plane  z  =  0  represents  the  free  surface.  Otherwise  the  foil  may  be 
placed  anywhere  along  the  z-axis. 

Function  XLE.  Function  XLE(N,S)  defines  the  x-coordinate  of  the  leading 
edge  of  the  foil.  Functions  N  and  S  are  the  section  number  and  fractional 
spanwise  coordinate  b' ,  respectively.  The  variable  NNNN,  placed  in  labeled 
common  block  CNRA,  provides  a  running  index  corresponding  to  succeeding 
configurations  treated  with  a  single  computer  run. 

Function  B.  Function  B(N,S)  defines  the  length  of  the  semichord  of  the 
foil  at  the  fractional  spanwise  position  S.  Function  N  is  the  number  of 
the  foil  section.  The  variable  NNNN,  placed  in  labeled  common  block  CNRA, 
provides  a  running  index  as  previously  described. 

Spanwise  Pressure  Modes 

End  Conditions.  As  described  above,  the  Widnall  program  predicts  load 
distributions  on  the  foil  by  summing  several  terms  of  a  series  which  is 
used  to  represent  the  pressure  jump  between  the  upper  and  lower  surface 
of  the  foil.  Each  term  of  the  series  consists  of  a  product  of  a  chordwise 
and  a  spanwise  pressure  mode.  The  chordwise  modes  are  taken  from  two- 
dimensional  wing  theory  and  are  built  into  the  program.  The  program  user 
must  specify  the  spanwise  pressure  modes. 

This  set  of  spanwise  pressure  modes  must  be  chosen  to  conform  to 
physical  flow  conditions.  In  particular,  each  set  of  modes  must  satisfy 
the  following  boundary  conditions  at  the  ends  of  foil  sections. 

a.  Foil  loading  drops  to  zero  at  the  end  of  a  cantilevered  foil  and 
at  the  free  surface  on  a  surface-piercing  foil. 

b.  Foil  loading  is  nonzero  at  an  intersection  of  two  or  more  lifting 
surfaces,  and  at  an  intersection  of  a  lifting  surface  and  a  reflecting 
plane. 

c.  Pressure  jumps  along  a  closed  path  around  an  intersection  of  two 
or  more  foil  sections  must  sum  to  zero. 


Note  that  the  above  conditions  place  restrictions  on  the  magnitude  of  each 
pressure  mode  at  the  ends  of  foils.  When  zero  loading  is  required  at  the 
end  of  a  foil,  each  of  the  spanwise  pressure  modes  must  become  zero.  How¬ 
ever,  when  nonzero  loading  contributions  are  required,  at  least  one  of  the 
modes  must  be  nonzero.  Each  mode  must  be  specified  on  all  foil  sections, 
although  the  shape  of  the  mode  does  not  have  to  be  the  same  on  all  the 
sections,  and  each  mode  must  satisfy  condition  (c)  independently  of  the 
other  modes. 

Certain  foil  configurations  suggest  boundary  conditions  on  the  span- 
wise  slopes  of  the  modes,  in  addition  to  the  above  conditions  on  the  magni¬ 
tudes  of  the  modes,  at  the  ends  of  the  foils: 

d.  Loading  on  submerged  foil  tips  is  approximately  elliptical. 

e.  Loading  near  the  free  surface  on  surface-piercing  foils  approached 
zero  less  steeply  than  does  an  elliptical  distribution. 

f.  The  load  distribution  slope  at  foil  intersections  is  approximately 
zero  (an  end-plate  effect). 

It  is  convenient  to  satisfy  all  of  the  above  boundary  conditions 
which  apply  to  a  given  foil  by  specifying  one  pressure  mode  having  all 
requisite  end-point  magnitudes  and  slopes.  Additional  modes  are  then 
formed  from  the  first  mode  multiplied  by  integral  powers  of  S  and/or 
(1.0-S).  The  additional  modes  permit  refinement  of  the  pressure  distribu¬ 
tion  by  the  calculation  process. 

The  pressure  modes  must  satisfy  one  final  condition: 

g.  Both  antisymmetric  and  symmetric  pressure  modes  must  be  included 
on  foil  surfaces  which  do  not  have  purely  symmetric  loading. 

Inclusion  of  both  symmetric  and  antisymmetric  modes  allows  the  program 
complete  flexibility  in  proportioning  the  loading  among  different  surfaces. 
A  commonly-encountered  foil  configuration  requiring  the  use  of  both  types 
of  modes  is  the  inverted-T  foil  at  a  nonzero  roll  angle  beneath  a  free 
surface,  an  example  of  which  is  treated  below. 
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Subroutine  FUNCTN.  The  subprogram  Subroutine  FUNCTN  (N,S,F)  is  used  to 
define  the  pressure  modes  for  use  in  the  program.  Function  F  is  an  array 
which  contains  the  ordinates  of  the  pressure  modes  as  functions  of  section 
number  N  and  fractional  spanwise  position  S. 

The  direction  of  a  pressure  jump  which  is  represented  by  a  given  mode 
is  designated  by  a  sign  convention  based  on  the  definition  of  the  local 
foil  slope  PSI.  The  positive  lift  direction  is  given  by  a  vector  n  normal 
to  the  foil  and  having  a  positive  z-component,  as  shown  in  Figure  B.3. 

This  convention  allows  interpretation  of  the  sign  of  the  lift  distribution 
output . 

Data  Card  Input 

Input  Format.  Data  cards  are  to  be  punched  with  the  input  quantities  shown 
in  Table  B.5  for  each  of  the  three  programs.  The  quantities  were  defined 
in  Table  B.4.  The  required  formats  are  shown  to  the  right  of  the  input 
quantities,  and  the  number  of  the  column  farthest  to  the  right  in  the  field 
in  which  the  input  quantity  is  to  be  placed  is  shown  above  the  quantity. 

Downwash  Input.  The  downwash  is  the  velocity  which  must  be  added  to  the 
free  stream  velocity  to  produce  the  flow  at  the  surface  of  the  foil.  The 
downwash  on  a  given  foil  surface  is  normal  to  the  foil  surface  and  to  the 
free  stream  (in  accordance  with  the  assumption  of  small  angles  of  attack) 
and  has  a  direction  opposite  to  the  direction  of  positive  lift  (n) .  Down- 
wash  magnitude  is  determined  by  the  boundary  condition  that  fluid  cannot 
penetrate  the  surface  of  the  foil.  In  the  case  of  steady  loading,  the 
boundary  condition  requires  tangential  flow  along  the  surface,  so  that  the 
downwash  to  be  input  to  the  program  is  v^/U,  where 

v 

—  =  sin  a  =  a  (Steady  Downwash)  (B.l) 
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TABLE  B. 5  -  DATA  CARD  INPUT 

TABLE  B . 5a  -  PROGRAM  1 
CARD  INPUT 


Program  1  (YKWS)  Sceady  Loading 


5  10  15 

NCASE/NCMBR/ IPUNCH/ 

Title  (any  80  characters) 

10  20  25  30  35  AO  A5  50  55  65 

RF/EM/NS/NOCP/NOLT/NOST/NCP/MP/ J3/ ARM/ 

5  10 

NI(2)/NI(3)/...NI(2*NS+2)/ 

10  20 

AF(1)/AF(2)/...AF(NS)/ 


5  15  25  35 

NOCP  Control  Point  Cards:  N/XOC/SOS/ETA 
Next  Card:  Downwash  Title  (any  80  characters) 

12  36  60 

(NOCP/ 3)  Downwash  Cards:  VR(1)/VR(2)/VR(3)/ 

. . . VR(NOCP)/ 

10 

Next  Ccrd(s):  YS/ (one  YS  card  for  each 
section — NS  total) 

(begin  at  Card  2  again  for  next  case) 


F0RMAT(3I5) 

FORMAT ( 1H1 , 80H) 

FORMAT(2F10. 5, 7I5.F10. 5) 

FORMAT (201 5) 

FORMAT  (10F10.5) 

FORMAT(15,3F10.5) 

FORMAT (1H1 , 80H) 

FORMAT (3(F12. 3, 12X) ) 

FORMAT(F10. 5) 


TABLE  B.5  (Continued) 

TABLE  B. 5b  -  PROGRAM 
2  CARD  INPUT 

Program  2  (YKWU)  Unsteady  Loading 


Card  1: 


Card  2: 

Card  3: 

Card  A: 

— 1—  Card  5: 


NCASE 


5  10  15 

NCASE/NCMBR/ IPUNCH/ 

Title  (any  80  characters) 

10  20  25  30  35  A0  A5  50  55  65 

RF/EM/NS/NOCP/NOLT/NOST/NCP/MP /.1 3/ ARM/ 

5  10 

N1(2)/NI(3)/...NI (2*NS+2)/ 

10  20 

AF(1)/AF(2)/...AF(NS)/ 


Sets 


5  15  25  35 

NOCP  Control  Point  Cards:  N/XOC/SOS/ETA 

Next  Card:  downwash  Title  (any  80  characters) 

12  2A  36 

(NOCP/ 3)  Downwash  Cards:  VR( 1 ) /VI ( 1 ) VR(2 ) / 

A8  60  72 

VI (2)/VR(3) /VI ( 3) / . . . VR(NOCP)/VI (NOCP) / 


(begin  at  Card  2  again  for  next  case) 


FORMAT (315) 
FORMAT(1H1,80H) 

FORMAT(2F10. 5,715, F10. 5) 

FORMAT (201 5) 

FORMAT ( 10F10.  5) 

FORMAT ( 1 5, 3F10. 5) 

FORMAT (1HI,80H) 

FORMAT(6Fl 2 . 3) 
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TABLE  B . 5  (Continued) 


TABLE  B.5c  -  PROGRAM 
3  CARD  INPUT 

Program  3  (YKWA)  Auxiliary  Program 


Output  I 
f  rom  . 
YKWS  \ 
or  I 
YKWU  I 


Card  1:  Title  (any  80  characters) 

Card  2:  Free  Surface  Card 

10  20  25  30  35  40  45  50  60 

Card  3:  RF/SPAN/NS/NOLT/NOST/NOCP/MP/NCP/EM/ 
65  70 

NCI-lBR/ARM/ 

10  20 

Card  4:  AF(1)/AF(2)/ . . .AF(NS)/ 

5  10 

Card  5:  N1 (2) /NI (3) / . . .NI (2*NS+2) / 


I  (PR]  matrix 

^  [PI]  matrix  (if  unsteady) 

l  Next  Card:  Downwash  Title  (any  80  characters) 


NCMBR 

Sets 


12  24  36 

(N0CP/3)  Downwash  Cards:  VR(1) /VI (1)/VR(2)/ 

48  60  72 

VI(2)/VR(3)/VI(3)/...VR(NOCP)/VI(NOCP)/ 
if  steady: 


10 

Next  Card(s):  YS/(one  YS  card  for  each 
section — NS  total) 


FORMAT (80H) 
FORMAT (80H) 


FORMAT(2F10. 5,615, 
F10. 5 , 15, F10. 5) 

FORMAT (6F10. 5) 
F0RMAT(14I 5) 


FORMAT (1 HI, 80H) 


F0RMAT(6F12. 3) 


FORMAT(F10. 5) 


and  a  is  the  local  angle  of  attack.  Note  that  downwash  velocity  is  given 
in  units  of  the  free  stream  velocity.  As  indicated  in  this  expression,  the 
downwash  input  to  the  program  is  the  local  angle  of  attack  (in  radians), 
in  keeping  with  the  theoretical  assumption  that  foil  angles  relative  to 
the  flow  are  small.  Camber  is  thus  represented  by  variations  in  vn/U  =  a 
as  a  function  of  x. 

Steady  downwash  is  input  to  program  YKWS  or  YKWA  using  the  VR  array. 
For  example,  a  downwash  input  of 

VR(I)  =  0.1,  1=1,  NOCP  (B. 2) 

would  indicate  an  uncambered  foil  at  an  incidence  of  0.1  radian. 
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Unsteady  loading  results  from  either  foil  oscillation  or  from  un¬ 
steady  flow  past  the  foil  such  as  that  due  to  the  orbital  velocities  of 
waves.  In  the  first  case  (foil  oscillation),  since  the  particles  follow 
the  motion  of  the  foil,  the  downwash  velocity  is  the  velocity  of  the  foil 
perpendicular  to  the  mean  position  of  the  foil,  which  is  assumed  to  be  at 
zero  angle  of  attack.  Downwash  input  is  the  sinusoidal  amplitude  of  the 
downwash  velocity,  in  units  of  the  free  stream  speed.  If  the  position  of 
the  foil  surface  is  given  as 

f (x,y , z, t)  =  0  (B. 3) 


then  the  normal  velocity  of  the  surface  is 


Dt  9t  9x  9t  9y  9t  9z  9t 


In  this  expression,  9y/9t  =  9z/9t  =  0  because  flow  is  restricted  to 
the  x-direction,  while  9x/9t  =  U,  the  flow  speed.  Since  foil  motion  is 
positive  upward  and  downwash  v^  is  positive  downward, 


=  _Df  =  _9f_  9f 

vn  Dt  9t  9x 


(B.  5) 


The  downwash  may  be  expressed  in  terms  of  the  amplitude  of  displacement  of 
the  foil  from  its  mean  position,  which  is  a  sinusoidally  varying  function 


f  =  f (x,y,z)  +  6(x,y,z,t) 


(B.  7) 
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Using  6,  the  downwash  becomes 


■  X  -LIT  °  iCJt 

1106  +U  -r —  e 
o  dx  I 


(B.  8) 


If  the  reduced  frequency  k  is  introduced,  the  downwash  may  be  expressed 
as  follows: 


k  =  coj/2U 


(B.  9) 


’  i2kU6  96 
o  ...  o 

v  =  -  -  +u  -5 — 

n  c  9x 


(B.10) 


which  is  then  normalized  in  terms  of  U  to  give 


6  36 

vn  .  __o  0  itot  Unsteady  Downwash 

U  1  c  3x  e  (Foil  Oscillation) 


(B. 11) 


The  real  component,  -96q/9x,  is  the  rate  of  change  of  the  unsteady  dis¬ 
placement  amplitude  in  the  streamwise  direction,  and  is,  therefore,  nonzero 
for  pitching  motion  and  zero  for  bending  motion. 

The  second  cause  of  unsteady  downwash  is  time-varying  flow  past  the 
foil.  An  analysis  is  presented  for  the  time-varying  flow  experienced  by  a 
foil  traveling  beneath  waves  due  to  a  head  sea  in  deep  water.  Surface 
waves  cause  the  water  particles  beneath  the  surface  to  move  in  circular 

paths  with  an  orbital  velocity  U' .  The  magnitude  of  the  orbital  velocity 
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at  a  given  depth  h  is  given  by 


U'  =  uj'  exp(2irh/A) 


(B. 12) 
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Both  vertical  and  horizontal  components  of  the  orbital  velocity  contribute 
to  the  unsteady  downwash,  but  at  a  relative  phase  difference  of  90  deg. 

The  unsteady  downwash  due  to  the  orbital  velocity's  vertical  component  is 


ulei(u>'t+TT/2) 


(B. 13) 


The  unsteady  downwash  due  to  the  orbital  velocity's  horizontal  component  is 


U'a  e 
o 


ico' t 


(B. 14) 


where  otQ  is  the  local  slope  of  the  foil,  whether  due  to  camber,  angle  of 
attack,  or  both.  The  total  downwash  due  to  wave  motion  is  then 


Vn  =  [u'a0+iu']  e1U),t  (B.  15) 

In  this  expression,  U'  and  aQ  are  to  be  evaluated  according  to  the  wave 
characteristics  and  the  location  on  the  foil  of  each  control  point.  The 
resultant  downwash  is  similar  in  form  to  the  downwash  for  pure  pitching 
oscillation  when  is  nonzero  and  is  similar  to  the  downwash  for  pure 
heaving  motion  when  a is  zero. 

The  computer  program  requires,  as  input,  values  for  the  unsteady  down- 
wash  seen  by  the  foil  traveling  through  the  waves  at  speed  U.  The  above 
expression  gives  the  amplitude  of  the  downwash  seen  by  the  foil,  but  u>' 
must  be  replaced  by  u),  the  wave  encounter  frequency.  The  resulting  down- 
wash  experienced  by  the  foil  is 


Vn  _  iwt  Unsteady  Downwash 

U  U  o  1  U  6  (Wave  Motion) 


(B. 16) 
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The  downwash  is  input  to  the  computer  program  in  the  form 


VR(I)  =  real  component 
VI (I)  =  imaginary  component 


1,  NOCP 


where  the  unsteady  exponential  e1Wt  is  assumed  by  the  program  and  is  not  a 
direct  input.  The  wave  encounter  frequency  to  is  input  as  part  of  the  re¬ 
duced  frequency  RF. 


Input  Restrictions.  Limitations  on  input  parameter  values,  based  on 
program  structure  or  dimensioning,  are  indicated  in  Table  B.4.  Suggested 
input  values  are  listed  in  Tables  2  and  6. 


PROGRAM  OUTPUT  DESCRIPTION 

All  three  programs  generate  printed  output  giving  the  loading  on  the 
subject  foil  systems  and  the  input  that  was  used  to  generate  that  loading. 
The  output  may  be  summarized  as  follows. 

First,  all  input  data  except  the  downwash  values  are  printed  out. 
Second,  Programs  1  and  2  print  all  matrices  involved  in  the  matrix  inver¬ 
sion,  including  the  resulting  unit  matrix.  The  size  of  the  off-diagonal 
elements  of  the  unit  matrix  indicates  the  quality  of  the  inversion.  Third, 
the  resulting  pressure  matrix  ([PR]+i[PI])  is  printed.  This  matrix  is  con¬ 
tained  in  the  input  for  Program  3.  Finally,  downwash  values,  pressure  mode 
coefficients,  load  coefficients,  and  center  of  pressure  locations  are 
printed  by  all  programs.  A  chordwise  pressure  jump  distribution  is  also 
given  by  Programs  1  and  3. 

The  pressure  mode  coefficients  appear  in  the  printed  output  as  the 
CR  and  Cl  arrays.  Array  CR  contains  the  steady  loading  coefficients  or 
the  real  part  of  the  unsteady  coefficient,  while  Cl  contains  the  imaginary 
part  of  the  unsteady  coefficients.  The  matrix  elements  correspond  to  the 
coefficients  a  found  in  Equation  (2)  in  the  following  manner: 
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Re  (a  )  =  CR(n+l,nri-l)  ' 

nm  n  =  0,1,2,... NOLT 

> 

m  =  0,1,2, . . .NOST 

Im(a  )  =  CI(n+l,m+l) 
nm  / 

Load  coefficients*  are  presented  as  both  spanwise  distributions  and 
total  loading  on  each  foil  section.  The  spanwise  distributions  consist  of 
lift,  pitching,  and  roll  moment  coefficients  per  unit  span  which  are  given 
at  intervals  of  5  percent  of  section  span  on  each  foil  section.  These 
spanwise  coefficients  are  defined  in  Table  B.6.  The  spanwise  pressure 
modes  (F)  selected  by  the  program  user  are  also  listed  at  5  percent  inter¬ 
vals  for  reference.  Total  lift  and  moment  coefficients  are  obtained  by 
numerical  integration  of  the  spanwise  distributions.  Unsteady  load  coef¬ 
ficients  are  given  in  terms  of  both  magnitude  and  phase  angle.  An  example 
of  printed  output  from  each  program  is  given  along  with  the  sample 
calculations. 

Punched  output  can  be  obtained  from  Programs  1  and  2  and  is  required 
for  running  Program  3.  When  a  value  of  1  is  input  for  IPUNCH,  all  neces¬ 
sary  parameters  for  running  Program  3  except  the  downwash  cases  are  punched 
by  having  them  written  on  the  punch  tape. 

A  center  or  pressure  is  calculated  for  each  section  whether  the 
loading  is  steady  or  unsteady.  For  steady  loading,  the  conventional  center 
of  pressure,  at  which  the  lift  appears  to  act,  is  calculated  and  is  given 
in  terms  of  its  x-coordinate  and  the  fraction  of  the  section  span  from  the 
end  of  the  section  at  which  S  =  0.  For  unsteady  loading,  the  center  of 
pressure  is  calculated  in  the  following  manner.  Let  the  lift  and  pitching 
moment  be  given  by 


.  COS  U)t 

0 

(B. 17) 

COS  (cot— 4>) 

(B. 18) 

*Load  coefficients  are  only  valid  on  sections  having  constant  slope 
(PSI  =  constant);  this  restriction  does  not  apply  to  pressure  jump  dis¬ 
tribution  coefficients  CR  and  Cl. 
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TABLE  B. 6  -  OUTPUT  LOAD  COEFFICIENTS 


CLR.CLI 

Real  and  imaginary  components  of  spanwise  lift  coefficient; 
(lift  on  a  unit  span)/qcQ;  positive  lift  opposite  to  \ 

positive  downwash 

CLMG 

2  2 

Magnitude  of  spanwise  lift  coefficient;  (CLR)  +  (CLI) 

CLPH 

Phase  angle  by  which  lift  leads  displacement;  tan  1 
(CLI/CLR) ,  deg 

CPMR.CPMI 

Real  and  imaginary  components  of  spanwise  pitching  moment 
coefficient  calculated  about  x  =  ARM,  positive  tending  to 
twist  foil  nose  in  direction  of  positive  lift;  (pitching 

moment  on  a  unit  span)/qcQ2 

CPMG 

Magnitude  of  spanwise  pitching  moment  coefficient; 

(CPMR) 2  +  (CPMI)2 

CPPH 

Phase  angle  by  which  pitching  moment  leads  displacement; 
tan"1  (CPMI/CPMR) ,  deg 

CRMR.CRMI 

Real  and  imaginary  components  of  spanwise  roll  moment  co¬ 
efficient,  calculated  about  section  root  chord  (where  S  =  0) ; 
positive  tending  to  roll  foil  in  direction  of  positive  lift; 

2 

(roll  moment  due  to  a  unit  spanl/qc^ 

CRMG 

2 

Magnitude  of  spanwise  roll  moment  coefficient;  (CRMR)  + 
(CRMI)2 

CRPH 

Phase  angle  by  which  roll  moment  leads  displacement; 
tan"1  (CRMI/CRMR) ,  deg 

CP 

Pressure  jump  coefficient  (equivalent  to  lift  per  unit  area 

2 

in  units  of  dynamic  pressure);  2Ap/pU 

where  the  moment  is  calculated  about  an  axis  at  x  =  a.  By  expanding  the 
moment  expression,  it  is  possible  to  express  the  moment  in  terms  of  com¬ 
ponents  which  are  in-phase  and  90  deg  out-of-phase  with  the  lift. 


M  =  M  cos  d>  cos  wt  +  M  sin  d)  sin  wt 
o  T  o 


(B. 19) 


1 
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The  distance  between  the  specified  moment  axis  at  x  =  a  and  the  unsteady 
center  of  pressure,  x  =  b,  is  given  by  the  moment  component  in  phase  with 
the  lift  divided  by  the  lift,  or 

b  -  a  =  (Mq  cos  <t>)/Lo  (B .  20) 

The  following  considerations  may  help  the  reader  to  understand  the 
nature  of  the  unsteady  center  of  pressure.  If  the  unsteady  pitching  moment 
were  calculated  about  an  axis  passing  through  x  =  b,  its  magnitude  would 
be  a  minimum.  In  addition,  it  would  have  no  component  in-phase  with  the 
lift  and  would,  therefore,  be  90  deg  out-of-phase  with  the  lift.  Measure¬ 
ment  of  either  of  these  characteristics  would  permit  a  direct  experimental 
determination  of  the  unsteady  center  of  pressure. 

The  program  output  gives  the  location  of  the  unsteady  center  of  pres¬ 
sure  in  terms  of  its  x-coordinate,  b,  and  the  fraction  of  the  section  span 
from  the  root  of  the  section. 

OPERATING  INSTRUCTIONS 

The  Widnall  lifting-surface  program  has  been  adapted  to  run  on  the 
Control  Data  Corporation  6700  digital  computer  at  DTNSRDC.  All  programs 
and  subprograms  are  written  in  the  FORTRAN  Extended  language  (Version  3.0). 

In  its  present  form,  the  program  requires  about  205,000  (octal)  memory 
locations.  This  amount  of  memory  permits  selection  of  the  full  range  of 
numerical  parameter  values  shown  in  Table  B.5.  The  standard  input  used  in 
the  present  study  could  be  run  with  substantially  less  memory  if  the  user 
reduced  the  appropriate  array  dimensions  in  the  program.  Conversely, 
larger  arrays  and  increased  dimensions  would  be  required  if  the  user  de¬ 
sired  more  control  points  and  larger  numbers  of  integration  intervals  than 
are  given  in  Table  B.4. 

Program  run  time  varies  according  to  the  value  of  numerical  input 
parameters.  With  the  standard  input,  an  execution  time  of  about  35  sec 
was  obtained,  expressed  in  terms  of  the  DTNSRDC  convention  of  time  required 
on  a  CDC-6400  computer.  Unexpectedly  and  for  unknown  reasons,  runs  made 


at  AR  =  1  required  about  twice  as  much  time  for  execution  (about  80  sec) 
as  runs  at  AR  =  6  and  10.  Compilation  of  the  source  deck  took  18  sec.  It 
is  concluded  that  the  computer  costs  for  using  the  program  will  normally  be 
small  compared  to  the  cost  of  learning  how  to  use  the  program  and  preparing 
input  data. 

SAMPLE  CALCULATIONS 

Sample  calculations  are  presented  to  illustrate  the  use  of  each  of 
the  three  programs.  All  data  and  subprogram  codings  corresponding  to  the 
sample  calculations  are  included  in  the  program  listings  given  in  Appendix 
C.  Discussion  of  each  case  consists  of  description  of  the  input  data  and 
presentation  of  printed  output  obtained  from  a  computer  run. 

Program  1  (YKWS) 

The  foil  configuration  chosen  as  an  example  for  the  steady  loading 
test  case  permits  several  aspects  of  program  use  to  be  illustrated.  How¬ 
ever,  the  calculated  loading  cannot  be  considered  accurate  because  the 
numerical  calculation  has  probably  not  achieved  a  stable  condition.  Con¬ 
sidering  that  60  control  points  were  recommended  on  a  single  foil  section, 
the  sample  calculation's  use  of  only  16  on  each  section  of  a  total  of  48 
on  all  three  sections  is  probably  inadequate  for  numerical  stability.  The 
stability  of  multisection  foil  calculations  has,  in  fact,  not  been  studied. 
Therefore,  this  sample  calculation  must  be  regarded  as  a  test  case  only, 
for  acquainting  the  user  with  program  capabilities  and  providing  a  means 
of  verifying  other  copies  of  the  program  at  relatively  low  cost. 

The  subject  foil  is  an  inverted-T  configuration  which  corresponds  to 
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a  hydrofoil  model  for  which  loading  has  been  determined  experimentally. 

The  foil  system,  shown  in  Figure  B.4,  is  rolled  relative  to  the  vertical 
and  is  surface-piercing.  As  a  result,  nonsymmetrical  loading  will  occur 
on  the  foils  even  though  the  strut  has  zero  angle  of  attack. 

Three  sections  are  used  to  represent  the  three  foil  surfaces.  Values 
of  NI(3),  NI(5),  and  NI(7)  have  been  set  equal  to  zero  to  prevent  inte¬ 
gration  over  negative  values  of  S  for  all  sections. 

Foil  geometry  is  specified  in  the  user-coded  subprograms.  Subroutine 
STRUCT  gives  the  y-  and  z-coordinates  for  each  section.  Leading-edge 
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Figure  B.4a  -  Foil  Planform 


4. 

35 

•-1.44-* 

Figure  B.4b  -  Strut  Planform 


Figure  B.4c  -  Foil  Configuration 
Looking  Upstream 


Figure  B.4  -  Foil  Configuration  for  Steady  Loading  Test  Case 


x-coordinates  are  given  in  Function  XLE.  Local  chord  lengths  for  the 
tapered  foils  and  the  untapered  strut  appear  in  Function  B. 

Spanwise  pressure  modes  are  specified  in  Subroutine  FUNCTN.  Ellipti¬ 
cal  loading  has  been  specified  at  the  tips  of  the  foils,  and  parabolic 
loading  specified  at  the  surface-piercing  end  of  the  strut.  Both  symmetric 
and  antisymmetric  modes  are  required  on  the  foils  in  order  to  permit  un¬ 
equal  loading.  Only  the  minimum  number  of  pressure  modes  has  been  speci¬ 
fied;  additional  modes  would  allow  refinement  of  the  spanwise  loading  dis¬ 
tribution  shapes. 

Control  points  were  placed  in  networks  of  4  spanwise  by  4  chordwise 
on  each  section,  making  a  total  of  48  control  points.  Downwash  values 
were  based  on  the  NACA  16-407,  a  =  1.0  foil  profiles,  with  the  uncambered 
strut  placed  at  zero  incidence. 

Printed  output  from  the  computer  run  is  shown  in  Appendix  D. 

Program  2  (YKWU) 

A  relatively  simple  rectangular  foil  configuration^  is  used  to  il¬ 
lustrate  an  unsteady-loading  calculation.  The  fully-wetted  foil,  shown  in 
Figure  B.5,  undergoes  forced  deflection  in  an  approximate  first  torsion 
mode. 

A  single  section  represents  one  half-span  of  the  foil,  while  S- 
symmetry  is  invoked  to  represent  the  other  half-span  produced  by  a  root 
reflection  plate.  The  symmetrical  half-span  is  included  in  the  integration 
by  making  NI(3)  nonzero.  Foil  geometry  is  specified  in  the  subprograms 
listed  in  Appendix  C.  Four  elliptical  spanwise  pressure  modes  were  se¬ 
lected.  Input  data  includes  the  reduced  frequency  value  of  1.97,  locations 
of  16  control  points,  and  unsteady  downwash  values  corresponding  to  the 
experimentally  determined  torsional  mode  shape  at  reduced  frequency  1.97. 
Again,  the  use  of  16  control  points  rather  than  the  60  recommended  for 
stability  is  intended  only  to  acquaint  the  user  with  the  program. 

Printed  output  is  shown  in  Appendix  E.  In  addition  to  the  printed 
output,  the  PR  and  PI  matrices,  along  with  several  auxiliary  quantities, 
were  placed  on  punch  cards  for  use  in  Program  3.  The  sample  calculation 
for  Program  3  uses  this  punched  output. 
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Figure  B.5  -  Foil  Configuration  for  Unsteady  Loading  Test  Case 


142 


Program  3  (YKWA) 

The  unsteady  load  calculation  described  above  is  continued  as  the 
sample  calculation  for  Program  3.  The  foil  geometry  and  pressure  mode 
subprograms  are  used  in  the  same  form  as  in  Program  2.  Data  cards  include 
the  punched  output  from  Program  2  and  the  same  unsteady  downwash  values 
used  for  Program  2. 

Printed  output  is  shown  in  Appendix  F. 
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Program  1  Listing  (continued) 
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Program  1  Listing  (continued) 
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SUBROUTINE  STRUCT (N.S.Z*Y.PSI ) 

C  HYSTAD  INVFRTFD  T-FOTL  STEADY  LOADING  TEST  CASE 
c  H=DFPTH  TO  FOIL  MIDSPAN  th=roll  angle  in  raoians 

C  NOTE  THAT  H  nrFINFD  POSTTIVF  AS  DEPTH/MN  CHD 
COMMON  /CNRA/  MNNN 
H=3.S? 

TH=0.174S1 
H=H* 1.301 
IF (N-2) 1  *?*3 
1  PSI=-TH 

Y=1 .R4Q5*S*C0S(TH)-m*TAM(TH) 

Z=-l ,949S*S*SIN(Th) -h 
RETURN 
?  PS I =-TH 

Y=-l ,P4RS*S*C0S (TH) -H*TAN (TH) 

Z=] .Q4RS*S*SIN(TH) -h 
RFTURN 

3  PST=(3.I41SR/?.n)-TH 
Y=-H*TAN(TH)*(1 .0-S) 

Z=-H*(l .0-S) 

RFTURN 

ENO 

FUNCTION  XLE(N.S) 

C  HYSTAD  INVFRTFD  T-FOTL  STFADY  LOADING  TEST  CASE 
COMMON  /CNRA/  MNNN 
IE (N-?) 1*1.3 

1  XLF=-1.0*  U .9405*0. 700?*0. 350) *AHS(S) 

RFTURN 

3  XLF=-0.4R3 
RFTURN 
END 

FUNCTION  P(N*S) 

C  HYSTAD  INVERTED  T-FOIL  STEADY  LOADING  TFST  CASE 
COMMON  /CNRA/  NNNN 
IF (N-?) 1.1.3 
1  B=1  •  0-0  •f>RS*S 
RFTURN 
3  b=1.44? 

RETURN 

FND 


Program  1  User-Coded  Subprograms 


SUBROUTINE  FUNCTN (N.S*F) 

e  hystao  invfrtfd  T-FOIL  STEADY  LOADING  TFST  CASE 
COMMON  /CNRA/NNNN 
OIMFNSION  F (6) 

DO  SIR  J-1,6 
SI 9  F(J)=n.O 
F=1 .0-S**? 

IF (N-2)  1*2*3 

1  F(1)=SOPT(F)*(-O.OS) 

F (?) =SOPT (F) * (♦n.OS) 

RFTURN 

2  F ( 1 ) =SQRT ( F ) *  (  ♦ 0  •  04  ) 

F(9)=S0RT(f )*(*0.0S) 

RFTURN 

3  F (1 )=F»0.09 
RFTURN 

Eun 


Program  1  User-Coded  Subprograms  (continued) 


1  1  1 

C  HYSTAD  INVERTED  T-FOIL  STEADY  LOADING  TEST  CASE-PROGRAM  CHECK  ONLY 
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Program  1  Punch  Card  Input 


3  .75  .9397  .0603 

3  .75  .7660  .1 

3  .75  .500  .1 

3  .75  .1730  .1 

3  .9698  .9397  .0603 

3  .9698  .7660  .1 

3  .9698  .5  .1 

3  .9698  .1730  .1 

C  HYSTAD  FOIL  AT  4  OFG.  INCIDENCE 

0.022  0.0?? 

0.022  0.062 

0.062  0.062 

0.096  0.096 

0.153  0.153 

0.153  0.02? 

0.022  0.022 

0.062  0.062 

0.096  0.096 

0.096  0.153 

0.153  0.153 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.000  0.000 

0.0 
0.0 
0.0 


10  10 
10  10 
10  10 
10  10 
10  10 
10  10 
10  10 
10  10 


0.022 

0.062 

0.096 

0.096 

0.153 

0.022 

0.062 

0.062 

0.096 

0.153 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


Program  1  Punch  Card  Input  (continued) 
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Program  2  Listing 
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Program  2  Listing  (continued) 
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SURROUTINF  STRUCT (N»S*Z»Y«PSI) 

C  UNSTEADY  LOADING  TEST  CASE— SWRI  RECTANGULAR  HYDROFOIL 
COMMON  /CNRA/  NNNN 
Y=S.0*S 
7=0.0 
PSI=0.0 
RETURN 
END 

FUNCTION  XLE(N.S) 

C  UNSTEADY  L0ADIN6  TEST  CASE— SWRI  RECTANGULAR  HYDROFOIL 
COMMON  /CNRA/  NNNN 
XLF=-1.0 
RETURN 
END 

FUNCTION  B (N * S) 

C  UNSTEADY  LOADING  TEST  CASE— SWRI  RECTANGULAR  HYDROFOIL 
COMMON  /CNRA/  NNNN 
R=1 .0 
RETURN 
END 

SURROUTINF  FUNCTN (N*S*F) 

C  UNSTEADY  LOADING  TEST  CASE— SWRI  RECTANGULAR  HYDROFOIL 
DIMENSION  F ( A ) 

R=SQRT ( 1 .0-S**?) 

F ( 1 ) =R 
F (?) =R*S**? 

F(3)=R*S*»4 

F(4)=R*S**ft 

RFTURN 

END 
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1  1  1 


C  UNSTEADY  LOADING 

tfst  case— 

SWt>T  RECTANGULAR 

HYDROFOIL 

1 

.97  0. 

0  1  16 

3  3  10 

10  0 

-0.38 

10 

10  10 

S.O 

1 

.20 

.?n 

.10 

1 

.40 

•  20 

.10 

1 

.60 

.20 

.10 

1 

.80 

.20 

.10 

1 

.20 

.40 

.10 

1 

.40 

.40 

.10 

1 

.60 

.40 

.10 

1 

.80 

.40 

.10 

1 

.20 

.60 

.10 

1 

.40 

.60 

.10 

1 

.60 

.60 

.10 

1 

.80 

.60 

.10 

1 

.20 

.80 

.10 

1 

.40 

.80 

.10 

1 

.60 

.80 

.10 

1 

.80 

.80 

.10 

C  SWRI 

rectangular 

HYDROFOIL  « 

TORSIONAL  MODE  ♦ 

RF=1 .97 

.176. 

-.035 

.176 

.104 

.176 

.243 

.176. 

.381 

.363 

-.072 

.363 

.215 

.363 

.501 

.363 

.787 

.559 

-.110 

.559 

.330 

.559 

.771 

.559 

1.211 

.770 

-.152 

.770 

.455 

.770 

1.062 

.770 

1.669 
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XCP>(SPHG/SLHG>*COS(  (SPPM-SLPnI*Q.0174S33IM-1./XZ>  ♦ARM 
YCP*(SRMG/SLMG)*CCS( (SRPH*SLRm>#0.01T4533>  MM./XZ) 

«Rm<6.ei<'i  xcp.ycp 

135  CONTINUE 

00  5u  1«1,N3CP 

sc  continue 

2  FORMAT (20H1CAM3ER  CASE  NUMBER ■ I S ( 3H  OF.JS.bn  CASES) 


ANStl  340 
ANStl 350 
ANSL1360 
ANStl 3  70 
ANSL1300 
ANSL1390 
ANSL1400 
ANStl4lO 


201  FORMAT  (//23H  CHECK  INPUT  PAPAMt  *  CRS/4X  ,  2MPF  ,  lx ,  2NEM  *  7X ,  2  2HNS  NCCP  ANSt 1420 


1N0LT  NOST  MCP,SX*2MHP,3X,3M6RM/  2F 10 . 5 ,€ 15* F 1 0 . 5  //7h  £  SPAN* 
2F1Q.SI 

203  FORMAT (1*0 ,39X*4HP£AL  *  11 X  ,4M  IHAG  ,  1 1 X  ,  wh  HAG  *  1  GX  ,5hPs*Sc./ 

13SH  SECTION  cIFT  COEFFICIENT  *  ,Fl0.fc,  2F IS . 6 , F IS. 2/ 

235M  SECT  ION  PITCH  NON  COEFF  ■  ,F10.E,  2F IS . b , F IS. 2/ 

335M  SECTION  ROLL  HOH  COEFF  «  ,F1C.«,  2F IS . 6 » F 15. 2> 

204  FORMAT (1H0.SX, 9*PR  MATRIX  //<4tl7.0»> 

2a5  FORMAT (1MQ,5X,9hPI  MATRIX  //(‘.Elf. 8)1 

20k  FORMAT (6F12. 3) 

230  FORMAT (l*0/(bFl2. 5)1 

2  05  FORMAT  <  1*  ,  J2X , 19HL0CAL  COEFFICIENTS 

ION  SECTION,  I2,23M  AT  SPANhISE  POSITION  S /// 

24bX,2 OHMOMENTS  CALCULATED  ABOUT  X  >  ,F6.2// 

3  34X*  5 IMPh ASE  IS  ANGLE  8Y  HHlC*  QUANTITY  tEACS  DISPLACEMENT 


ANSL 1 A  30 
ANSL1440 
ANSL14S0 
ANSL1460 
ANSL1470 
*NSH4»0 
ANSL1490 
ANSt 1500 
ANSL1S10 
A*.  St  1520 
1 CH ANStl 5  30 
ANSL15N0 
ANSL15S0 
IANSL1S60 


210  FORHAT  (1hO*2X* lHS*GX*3HCLR*7X>3HCLI,FX*AHCLMG»kX,4HCL*H,EX,AHCPPK,ANSLlS 70 

l6X,4MCPHI,bX,4MCPNG,6X,4MCPP*,bX,4MCRHk,6X,4NCRHI,6X,.»MCRN&t6X  ,4HCANSL1S00 
2RPH  /3X,lH-,6  ANSL1S90 

3X*3M *  7  X  *  3  H ,  7  X  ,  4H - -,b*  ,4H - *bX,4H - ,  6  X , 4H ,  bX,  AH - AKSL16  00 

4,bX,4M - -  ,bX,4M - ,bXt4H- - ,bX,4M - *SX,4H - )  ANSI  lb  10 

211  FORHAT  ( 1*  VFS.2*3(3F10.6*F9.2.1X> I  ANSL1620 

212  FORHAT  HMD  ANSL1630 

214  FORHAT <1MO,34HX-COOROINATE  OF  SECTION  CP  **  FlO.6/  ANSL1C 40 

13SH  SPAN  WISE  COORDINATE  OF  SECWCN  CP>  .FiO.fcl  ANSL1ES0 

215  FORMAT  <//// 3X , 1HS  *  SX *  ANSL1660 

2  4Hf  <11  ,5X,4MF<2> fSX,4NFI3) * 5X ,4FF  (4) *SX*4HF (5) ,SX  ,4HFANSL167| 

3(61  /  ANSL1M0 

4  31  »lH*t  SX,4M--  —  ANSL1690 

>,$X*4H—.|$XI<M - *5X,4M - ,5X,4M - ,SX,  4H - 1  ANStl  7  QQ 

21 1  FORMA T (1H  ,FS.2* 1P6C9.1)  ANStlflO 

219  FORMAT (1H  *  AH  SF(1I«  10F10.SI  ANSL1720 

220  FORHAT ( 1H  *OH  Nl(2)»  20IS)  ANSL173I 

230  FORMAT (1H  tlPSEl2.3l  ANSL1740 

2  31  FORMA  T(1H  ,///99n  REAL  COtfflClEtl  MAT***  CM  <SP*M»lSZ  M  CHO*DHSl  ANStl  750 


111 

241  FORMAT ( 1M  ,///S0H  IHAG  COEFFICIENT  MATRIX  Cl  (SPANMISE  X 
II) 

900  FORHAT  (1H|)*20A4) 

923  FORMAT (2Fi0.S*bI5,Fi0.S,IS,Fl0.5l 

924  FORMAT (4E17.AI 

925  FORMAT (bFi&. 5) 

92t  FORMAT (14151 

END 


ANStl 760 
CHORObISE ANSL1770 
ANSU1700 
ANSL1790 
ANStl A 10 
ANSL1010 
ANSt II 20 
ANSL1030 
ANSL1040 


Program  3  Listing  (continued) 


FUNCTION  XLE(N.S) 

C  UNSTEADY  LOADING  TEST  CASE — SWRI  RECTANGULAR  HYDROEOIL — RE=1.97 
XLE=-1 .0 
RETURN 
END 

FUNCTION  R (N»S) 

C  UNSTEADY  LOADING  TEST  CASE— SWRI  RECTANGULAR  HYDROEOIL — RE=1.97 
B=1 .0 
RETURN 
END 

SIJRROUTINE  EUNCTN  (N,S«E) 

C  UNSTEADY  LOADING  TEST  CASE  — SWRI  RECTANGULAR  HYDROEOIL— RE=I  .97 
DIMENSION  E (6 ) 

R=SQRT(1.0-S**?) 

E ( 1 ) =R 
F (?) =R*S**? 

E (3) =R*S**4 
E(A)=R*S**6 
RETURN 
END 
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r  UNSTEADY  LOADING  TFST  CASE — SWRI 
J3=0  INFINITE  FLUID 
1.97000  5. 00000  1  3  3 

1.00000 

10  10  10 

.1R477725F-01  -.30675943F-01 

.171F1?H3R*00  .6386P057F ♦OO 

-.752?7897F*00  -. 1 0 1 96953F-0? 

.96467961F-0?  .2447867PF-01 

-.P7941997E-01  .68 1 35495F*00 

.107F3080F>00 

-.345»0948F-01  .556SPP44F-01 

.2P146141F-01  -.10FF48?»F*00 

.58583096E*00  .23P1PP20F-01 

.33613705F-0?  - .49O7S540F-0? 

-.70764595E-01  -.2668932?E*00 

•  50662731 F  *00  .251 383?SF-03 

-.599S10«1F-0P  -.67P79098F-0? 

.1S9P4178F-01  -.4S8C(6014F*n0 

•171 35733E-0 1  - . 1 3SSS738F -0 1 

.P4R99114F-01  -.37P79304F-01 

.60652896F-0?  - . 1 76396P4F-Q 1 

-.P409P3PSF*00  -.156?69«0F-01 

-•6RRP27R5F-0?  .l?10q7Q7F-01 

-• 1909219PF *00  -.70P??94PF*00 

•9611 1445E*00  •  . 70063470F-03 

-.969P59P4E-0?  -.?61 1 GF?9F-01 

.368PP909F-01  -.P4867990E*00 

-.415163SSF-01  -.3P307PPOF-01 

.441698P7F-01  -.75778312F-01 

-.59804207E-02  . 3981 08S6F-0 1 

-.63497004E*00  -.27752678E-01 

•P473154RF-0?  -.4160SR70F-02 

.79S13627F-01  . 30P051 30F *00 

-.RP44R680E*00  -.PP54391 6F-03 

.7P596514F-0?  .26710346F-0? 

-.212P6610E-01  .7P17P1 1 9F*00 

.14168333F-01  . 149S6353F-01 

-.385P9186F-01  .6?^?3197F-01 

•25476R70F-0?  -. 1 3S6RR81F-01 

•251 96399E*00  .21 6056P6F-01 

•2784P759F-0?  .586490O9F-0? 


rectangular  hydrofoil 

16  10  10  0.00000  1 


-.11?28947E*00  -. 1 007 1 209E*00 

.1 1804733E*00  -.20441520E*00 

-.57613615F-02  . 10391 180E*00 

-.5776]464F*00  -. 1201 1244E-01 

-.45963P89E-02  .81337399F-02 

-.47P14427F-01  -.60335955E*00 

.71269003E*00  . 39861609E-02 

-.1994P488F-01  -• 1 00 16424 E ♦OO 

•  1 1458083E*00  -.69026990E*00 

-.17444346F-01  .39689233F-01 

-• 77143908E-0 1  . 1 3569095E *00 

-.16810617E-02  .17340131E-01 

.23824915E*00  . 1 0 1 P9240E-0 1 

-.40216593F-03  . 1 2 1 P0437E-02 

.19814147E-01  .25341759E*00 

-.485?0604E*00  . 83796944E-03 

•78610360E-02  . 30 1 0RR80E-0 1 

-.66090388E-01  .463R0663E*00 

.44706184E-01  . 1 1026120E*00 

-• 15152783E*00  .26216933E*00 

• 14143330E-02  -. 38674223E-0 1 

.61912887E*00  . 1 2673153E-0 1 

.21459717E-02  -.35967466E-02 

.54825850E-01  .65690008F*00 

-.90181716E+00  -.119R9277E-0? 

.P0P691 19E-01  . 1 0679774E*00 

-. 14874322E*00  .87267536E*00 

-.16094111F-01  -• 43653 185E -01 

. 12571R59E*00  -.22177948F*00 

-.39P52748E-03  . 1 3865961E-01 

- . P536 1 066F ♦OO  -. 15304997E-01 

-.65190029E-03  . 1 1991497E-02 

-.23406082E-01  -.26544781E *00 

.75P3R1 18F*00  .44277223E-03 

-.62931 141E-02  -.27661 187E-01 

.  1  0733603F.  ♦OO  -.  7257588  IE  *00 

.63458329E-01  -. 168451 16E-01 


Program  3  Punch  Card  Input 


I 


-.225637826-01 
.409456016*00 
.872998956-02 
•  1 34041  OOF *00 
.874752956-02 
-.353865276-01 
•55270247F-02 
• 34653487F ♦OO 
.357678086-03 
.338842986-02 
-.271817766*00 
-.44586R90F-02 
-.810371006-01 
-.108413016-02 
.254580876-01 
.236776706-02 
”•111 ?41 30F  *00 
“•1251 18916-02 
.250 1 240 1F-0 1 
-.513501576*00 
-.890902946-02 
-.174623866*00 
-.196268006-02 
.422058476-01 
-.107604986-02 


-.34785077F+00 
-. 17649Q03F-02 
1 1 60406 1 F ♦OO 
”  •  1  8®  255 1 6F-0 1 
.2972627OF-01 
. 14P99775F+00 

-.635*i664HF-01 
,?406P757F-01 
•21384147F-02 
•  140 1 1 575F  +  00 
-.P89S9528F-04 
.38866654F-OI 
-.50173080F-Q3 
-.135220646-01 
-.92297972F-01 
-.1381 P1P3F-01 
-.1498P27SF-01 
-.1PQ47438F-0? 
.373492806*00 
.58241 137F-03 
. 125P3551F+00 
.375 1^1 79F-0 1 
-.3O9095P4F-01 
-.19S12794F.Q0 
. 27420  0  33F-0 1 


. 19965080F-01 
.239000P2F-01 
. 1810«455F-01 
-.S2204462F.-02 
-. 12P33P76F*00 
•  75389699 F -01 
-.20345209F-01 
.201 1««28F-01 
•111 336026-01 
“•1351 7753F-0 1 
.58314003F-02 
.274670 66F-02 
71530756F-03 
.45507849F-01 
-.782166406-01 
.8087281 1F-02 
-.54109258F-02 
-.23170867F-01 
-.26128349F-01 
-.7419111 3F-02 
-.26731414F-01 
.194365276-02 
. 140871 066*00 
-.41532819F-01 
.256277226-01 


-.426597566*00  -. 355 1 3202F-0 1 


“ . 31 330  309F-0 1 


399359266-03 

2883021 3F-02 
42838977F  *00 
43514301F-02 
132273606*00 
1 1791485E-02 


571 27361F-03 
14690516F ♦OO 
18732584F-03 
38912777F-01 
29624449F-01 
12888224F-01 


8 1 346054E-02 
22635650F-01 
28652938F-02 
99 1 95763F-02 
68598156F-03 
478434 1 9E-0 1 


*  • 

36321571F-01 

.151882196*00 

.61592235F-01 

• 

47716782F-03 

-.944990206-02 

-.123786426-01 

• 

1 7250962F ♦OO 

.304932666-01 

.154645936-01 

C  SWRT  RECTANGULAR 

hyorofoil*  torsional  mqdf«  rf=i. 

.176 

-.035 

.176 

.104  .176 

.176 

.381 

.363 

-.072  .363 

.363 

.501 

.363 

.787  .559 

.559 

-.330 

.559 

.771  .559 

.770 

-.152 

.770 

.455  .770 

.770 

1.669 

97 


• 247939B0E-0 1 
.42646661E-02 
.10594658E-01 
•2531 4149E-01 
.5941 1 380E-01 
•2976781 OE-OP 
.19035698F-01 
. 39933747F ♦OO 
.72088548F-02 
. 1 0458343E-01 
.203421 07F -02 
.H1089504F-02 
. 5677258 0E-02 
.45399362F-01 
•29677536F-03 
.38955707F-03 
•?23299B2E*00 

•  1 9056987F-0 1 

•  3548940 1E-0 1 
.158489756-02 
.121580996-01 
• H54 1 3046F-02 
.347976606-01 
.177808296-02 
•21 748995E-0 1 

•  593369406'  ♦OO 
.861 16472E-02 
•20833747F-0 1 
• 73796541E-03 
. 1 1 300588F-0 1 
.314652666-02 
.24397777F-01 
.43786899E-03 
.13914201E-03 
.4749471 1E*00 

.243 

.215 

-.110 

1.211 

1.062 
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APPENDIX  D 

PROGRAM  1  OUTPUT  FOR  SAMPLE  CALCULATION 


. 172606096-01 
. 239000606-02 

.eoi^qmze-qj 

.6609  2966E -01 
.09372066E-Q3 
•  20296196E  -0  J 

16  96  6263E -0 1 
. 37992693E-02 
.6156<092E-03 
•  •  6  0  70  63 36E - 0  1 


.  17000619E-03 
•96966009E-C2 
.36623167E-Q2 
.201177UE-03 
.  12672906E-02 
-.71679779E-03 
•  255  3#  #Cfl£-0  3 
•60366621E-02 
•  »  S6910200E-02 
.396390626-03 


.3l($109lE-0J 
•  3065  129*E-01 
.226001906-02 
.616636296-03 
-.906016296-02 
.237936906-02 
.609069906-03 
-.36621270E-01 
.610662006-02 
.666267666-03 


•  166167016-02 
. 22  1660 95E- 03 
. 6  19760866-02 
.65012669E-O2 
• 1369220  7L-03 
. 366351 76E-02 
- . 1 960232  6E- 02 
.369697066-03 
,66  7372  066-02 
-.60993000E-02 


CAMBER  CASE  NOH6E9-  1  OF  1  CASES 
C  HTSTAO  FOIL  AT  6  OEG.  INCIDENCE 


.62200  .02200  .02200 

.02200  .06200  .06260 

•06200  .06200  .09600 

•09606  .69600  .09600 

.19300  .19360  .19300 

.19300  .02200  .02200 

.02200  .  02200  .  06260 

.06200  .06200  .06200 

.09600  .09600  .09600 

.09600  .19300  .19300 

.19300  .19300  0.00000 

0.0C000  0.00000  6.00000 

0.00000  0.60000  0.00000 

0.00600  0.00000  0.00000 

0.00000  0.00000  0.00006 

O.OCCOO  0.00000  0.  00000 


NEAL  COEFFICIENT  MATRIX  CR  (SPANNISE  X  CMCBONISE) 
•1.7611-09  1.732E-07  -2.160E-06 

1.21CE-02  1.213E-02  -1. 9966-02 


LOCAL  COEFFICIENTS  ANC  SPANNISE  FUNCTION  VALUES  ON  SECTION  1  AT  SPANNISE  POSITION  S 
MOMENTS  CALCULATED  ABOUT  X  >  -0.060 


S 

CLR 

CPMR 

CRNR 

F  (11 

F  ( 21 

- 

— - 

mmmm 

.... 

mmmm 

0.00 

.306702 

.060600 

0.000000 

-9.0E-02 

9.0C-02 

.09 

.306390 

.031326 

.016933 

-9.0E-02 

9.0E-02 

.10 

.309266 

.022021 

.029796 

-9.0E-02 

.19 

.  303311 

.012792 

.066360 

-6.9E-02 

.20 

.  J00906 

.003909 

.090999 

-6.9E-J2 

6.9E-02 

.29 

.297060 

-.009396 

•072309 

-6.OE-02 

6 .06-02 

•  30 

.292691 

-.01612J 

.009979 

-6.0E-02 

6.0E-02 

.39 

.207370 

-.022910 

.090063 

-6.7E-02 

.60 

.201170 

-.030699 

.109620 

-6.6E-02 

6.6E-02 

.69 

. 27  J969 

-.037960 

.120171 

-6.9E-02 

6.9E-02 

.90 

.269601 

-.066009 

.129606 

-6.3E-02 

.99 

.296  2  1  3 

-.090926 

.137399 

-6.2E-02 

.60 

.269629 

-.096167 

.163937 

-6.0E-02 

.69 

.233136 

-.060372 

.167711 

-3.0E-02 

3.0E-02 

.70 

.219006 

-.063320 

•169600 

-3.6E-02 

3.6E-02 

.79 

.202917 

-.066762 

.160369 

-3.3E-02 

3.3E-02 

.00 

.106069 

-.066207 

.163937 

-3.0E-02 

.09 

.161607 

-.061307 

.133090 

-2.6E-02 

2.6E-02 

.90 

.  13372! 

-.096796 

.117312 

-2.2E-P2 

2.2E-02 

.99 

.099793 

-.062106 

.000709 

-1.6E-02 

1.6E-02 

1.00 

.060000 

-.000000 

.000000 

-1.3E-00 

SECTION 

LIFT  COEFFICIENT 

■ 

.260933 

SECTION 

PITCH  HCH  COEFF 

■ 

-.029660 

SECTION 

ROLL  MCM  COEFF 

■ 

.099666 

X-COOPOINATE  OF  SECTION  CP  ■  .266276 

SPANNISE  COORDINATE  OF  SECTION  CP*  .027339 


Program  1  Output  --  Steady  Loading  Test  Case  (continued) 
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LOCAL  COEFFICIENTS  ANQ  SPANNISE  FUNCTION  VALUES  ON  SECTION  2  AT  SPANNISE  POSIT  ION  S 
NONCNTS  CALCULATED  A|OUT  I  •  0.000 


s 

a* 

CPNP 

CPPP 

Fill 

F 1 21 

... 

-•-» 

mmmm 

*••• 

0.00 

.  306290 

.0*0*6* 

Q.OOOOuC 

*.  0E-0  2 

5.0E-02 

.09 

.30906* 

.011209 

.01*907 

*  •  OE-O  2 

5 .06-02 

.10 

.  30*71* 

.021922 

.029702 

*. 0E-02 

5.0E-02 

.19 

.  302709 

.012672 

.0**271 

*. 0E-02 

* .9E-02 

.20 

.300062 

.003920 

.050*97 

3. 96-02 

* .9E- 02 

.29 

.296929 

-.005*17 

.072259 

3.9E-02 

* • 0E- 02 

.10 

.2921*3 

-.01*1*0 

.005*30 

3.OE-02 

* . OE*  02 

.19 

.206079 

-.022526 

.097072 

3.7E-02 

* • 7E- 02 

.*0 

.200602 

-.030*06 

.109*30 

3.7E-02 

*  «6E- 02 

.*9 

.271*90 

-.03793* 

.119963 

3.6E-02 

* .5E- 02 

.90 

.269220 

-.0**760 

.129262 

3.5E-02 

* . 3E- 02 

.99 

.299769 

-.050070 

.137121 

3.3E-02 

* .2E-02 

.00 

.2*9000 

-.056100 

.1*1200 

3.2E-02 

* . OE- 02 

.09 

.212729 

-.06029* 

.1*7*5* 

3.0E-02 

3 .OE- 02 

•  TO 

.210700 

-.0632*2 

.1*9229 

2.9E-02 

3.6E-02 

.79 

.202569 

-.06*671 

.1*0000 

2.6E-02 

1.3E-02 

.00 

.103790 

-.06*193 

.1*3200 

2.  *E-02 

3.0E-02 

.09 

.161127 

-.061215 

.133665 

2.1E-02 

2. EE-02 

.90 

.  13J*91 

-.05*669 

.117100 

1.7E-02 

2.2E-02 

.99 

.099626 

-.0*20*0 

.000551 

1.2E-02 

1.6E-02 

1.00 

.000900 

.OOOOuO 

1.0E-00 

1.3E-00 

SECTION 

LIFT  COEFFICIENT 

■ 

.2*0515 

SECTION 

PITCH  NON  COEFF 

• 

-.029652 

SECTION 

POLL  NCN  COEFF 

■ 

.099*93 

I-COOPOINATE  OF  SECTION  CP  •  .2*697* 

SPANNISE  COONOINATE  OF  SECTION  CP*  .027135 


LOCAL  COEFFICIENTS  AKO  SPANNISt  FUNCTION  VALUES  CN  SECTION  3  AT  SPANNISE  POSITION  S 
MOMENTS  CALCULATED  AeoUT  I  •  3.060 


S 

CLP 

CPNP 

CPNP 

Fill 

F 

• 

— - 

•— 

— — -» 

•  ••  • 

0.00 

-.000169 

.006027 

0.000000 

9.  OE  - 02 

0. 

.05 

-.  000160 

.000027 

-.000330 

9.0E-02 

0. 

•  10 

-.000165 

•000627 

-.000059 

0. 9E-02 

0. 

.15 

-.000161 

. 0  Oo  02  7 

-.000067 

O.OE-02 

0. 

.20 

-.00035* 

.000026 

-.00011* 

0.6E-02 

0. 

.29 

-.0001*6 

.006026 

-.000160 

0.*E-02 

0. 

.10 

-.000116 

.000025 

-.000163 

0.2E-02 

0. 

.15 

-.00012* 

.00002* 

-.000103 

0. 

•  *0 

-.000110 

.000023 

-. 00u200 

7 ! 6E -02 

0. 

.*5 

-.00029* 

.006022 

-.00021* 

7.2E-02 

0. 

•  50 

-.000277 

•060020 

-.006223 

6.7E-02 

0. 

.55 

-.000257 

•606019 

-.000220 

6. JE-02 

0. 

.60 

-.  000216 

•000017 

-.060229 

5.0E-02 

0. 

.69 

-.000213 

.000016 

-.000223 

5.2E-02 

0. 

.70 

-.000106 

.60001* 

-.000213 

*• 6E-02 

0. 

.75 

-. 000161 

.000012 

-.000195 

3.9C-02 

0. 

.00 

-.0001  J3 

.000010 

-.000171 

1.2E-02 

0. 

.05 

-.000102 

.000600 

-.0001*0 

2.5E-02 

0. 

.90 

-.  000070 

.000009 

-.000102 

1.7E-02 

0. 

.95 

-. 000016 

.000003 

-.000055 

O.OE-03 

0. 

1.00 

-. 000000 

.000000 

-.000000 

5.0E-K5 

0. 

SECTION  LIFT  COEFFICIENT  •  -.0002*6 
SECTION  PITCH  PCN  CCEFF  •  .000010 
SECTION  PCLL  NCP  CCEFF  ■  -.0001** 

l-COOPO IN AT E  OF  SECTION  CP  ■  .211*20 
SPANNISE  COONOINATE  OF  SECTION  CP*  1.7**129 


Program  1  Output  --  Steady  Loading  Test  Case  (continued) 
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Program  1  Output  —  Steady  Loading  Test  Case  (continued) 
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-.51350157£»C0 
-.81-9*2946-42 
• • 12482  386E ♦ J  3 
- •  1 162  68  0 06 -  02 
.422358476-31 
- • 1 3  75*  498E -02 
-.42S597‘»b£  *  JO 
.399359256-03 
-.259332l3E-b2 
.42338977E**., 
.435l43.'lF-)2 
. 1 32  2  736  0E ♦ C  0 
.117914856-C2 
-.3bJ21S71t-dt 
.47716782F-03 
.17250  962E*..- 


.  5  8241 13  7E  -  3  3 
.I2563551fc>00 
.375131796-31 
-•  309U9554E-01 
I95l279*e*0o 
•  2  7426  0  33E-0 1 
• «  355 13  2U2E-0 1 
• 57l273blt-03 
-.l45905lbE*00 
-.187325846-53 
- •  3  8 11277  7£ -0 1 
•  •  295244496-01 
. I2888224t-0l 
.  1S188219i>Ow 
-. 94499w2ut-02 
.3U493266E-Q1 


741911136-02 
-•  257314146-*  1 
. 19435527t-42 
.  14j871*5LU' j 
-.415328196-01 
.25527722E-CI 
-.3133j319E-*1 

•  61 J45054E-U2 
.22535554E-11 
. 286S2938E-Q2 

•  99195753E-4  2 
-.695981566-01 
-.4784  3419E-01 

. bl 592235E -0 1 
- • 12378642E-0 1 
. 154645936-01 


.  158481756-02 

•  1 21583  99E -0 1 
854130466-02 

.347976646-01 

-.177808296-02 

•  21 7489956-3 1 
.5933b940£»J0 

-.8bllb472£-02 
.248337476-01 
-.73796541E-03 
-.113405886-31 
. 31465266E-02 
-.243977776-31 

•  437868996-0  3 
.  13914201E-J3 

47494711E»U0 


COMBER  COSE  NUMBER*  1  OF  1  COSES 

C  SMRX  RECTANGULAR  MVOROFOIL,  TORSIONAL  MQOE,  RF*1.97 


.17500 

- 

.-»3 54J 

.17535 

.13443 

.17543 

.24300 

. 1760U 

. 38100 

.35300 

-.07200 

.36300 

.21500 

.  16  30  n 

.5414b 

.36340 

.78700 

•  5590*1 

-.11303 

•  5593  0 

.33040 

•  5594  C 

•  7710b 

.5590b 

1.21100 

.77,,ua 

- 

•  1 521 3 

.773*0 

•  455b  0 

. 7700- 

1 .36203 

•  77  3  m  j 

1 

.56913 

REOL  COEFFICI 

ENT 

MATRIX  CR 

(SPANMISE  x 

CHOROHlSE) 

4  •  Q90E- 03 

-1. 

0786-02 

1.855F-02 

2. 6286-32 

-9. 

5476-02 

1.88C6-31 

-3.373E-03 

9. 

373E-05 

3.3536-02 

IMOG  COEFFICIENT 

MATRIX  Cl 

<SPANWISE  X 

CMORObISEI 

1.703E-C4 

2. 

3  ■)  36  -  «  2 

1. 7516-04 

-2.C9*E-Q3 

1. 

7436-31 

2.2*46-33 

-4.981C-CJ 

-0. 

5*  7E-C3 

f.055£-fl3 

S 

F  1 11 

F  ( 2) 

F|3> 

0.00 

1. JE*30 

b. 

0. 

.05 

1 .  CE  +  iib 

2. 56-43 

6.2E-36 

•  16 

9.9E-01 

9.9E-03 

9.96-35 

.15 

9.96-ul 

2.2E-42 

5. 36-*  4 

.20 

9.8E-01 

3.96-42 

1. 66-03 

.25 

9.7E-01 

6.1E-C2 

3.8E-33 

•  30 

9.5E-C1 

8.6fc-«2 

7.76-vJ 

.35 

9.4E-01 

1.  16-01 

1.4E-U2 

•  40 

9.2E-01 

1.5E-01 

2. 36-02 

.45 

8.9E-01 

1.86-41 

3.7E-02 

.50 

8.76-01 

2.2E-01 

5.4E-02 

.55 

8.4E-01 

2.  5t-(i  1 

7.6E-02 

.60 

8. Jc-Ql 

2.9E-01 

1.4c- Jl 

.65 

7.6E-01 

3.2E-C1 

1.4E-H 

.76 

7.1E-01 

3.5E-C1 

1.76-01 

.75 

6.6E-01 

3.76-41 

2. 16-31 

.It 

6.0E-01 

3.8E-81 

2.56-01 

.85 

5.3E-01 

3. 86-01 

2.7E-41 

•  90 

4.4E-C1 

3.56-01 

2. 96-01 

.95 

3. lE-wl 

2.86-tt 

2.56-31 

l.ll 

2.5E-07 

2.56-07 

2.5E-07 

Program  2  Output  --  Steady  Loading  Test  Case  (continued) 
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LOCAL  COEFFICIENTS  ON  SECTION  1  AT  SPANNISC  POSITION  S 


s 

CLR 

CLI 

MOMENTS  CALCULATED  AMOUT  X  >  -.36 

PHASE  IS  ANGLE  9Y  WHICH  QUANTITY  LEAOS  DISPLACEMENT 

CLMG  clph  cpmr  CPNI  cpmg  CPPH 

CRNR 

CRMI 

CRNG 

CRPH 

- 

--- 

-- - 

.... 

.... 

— -- 

- -- - 

--«■■» 

— — 

---- 

—  — — 

0.00 

• . 25oS66 

2.3073  76 

2. 321332 

99.35 

.365066 

-.626659 

.563212 

-69.59 

.00k  5  00 

0.000000 

0.000000 

0.00 

.09 

>.297160 

2.366106 

2.361266 

99.50 

.372773 

-.636512 

.576023 

-69. 5g 

-. 333392 

.293263 

.295156 

96.56 

.&• 

-.297791 

2.962612 

2.660356 

96.90 

.395721 

-.659332 

.606266 

-69.25 

-.076666 

.615603 

.620066 

96.90 

•  19 

-.366699 

2.653590 

2.676379 

97.66 

.633566 

-.696695 

.659676 

-68.89 

-.  13w666 

.995096 

1.003661 

97.66 

•  20 

-.616966 

2.915616 

2.965523 

96.17 

.665761 

-.566692 

.732666 

-68.67 

-.209332 

1 .657809 

1.672762 

96.17 

•  29 

- • 50  7635 

3. 262  61 J 

3.262272 

98.69 

.55167- 

-.613113 

.626637 

-66. C3 

-.3171*,/ 

2.026756 

2.051628 

96.89 

.30 

-.913936 

3.627753 

3.679265 

99. 6U 

.629536 

-.669633 

.933611 

-67.6- 

-.<•60227 

2.720615 

2.759666 

99.60 

•  J9 

-.735726 

6.L61211 

6.127315 

1*3.27 

.716563 

-.775765 

1.  1 57626 

-67.19 

-.663762 

3.553559 

3.611601 

100.27 

•  AO 

-.671731 

6.57199 7 

6.615076 

1C3.69 

.616710 

-.670130 

1.193376 

-66.61 

-.671731 

6.531997 

6.615*76 

1«*0  .69 

.65 

-1.319121 

5.026760 

5.129067 

1J1.69 

.921627 

-.97JC10 

1.336166 

-66.66 

-1.166511 

5.655127 

6.770176 

101.66 

•  90 

-1.176997 

5. 52  96C  8 

5.653232 

101.99 

1.031179 

-1.072677 

1. 667796 

-66.12 

-1.666372 

6.912260 

7.066502 

101.99 

.99 

-1.336366 

6.022500 

6.196557 

102.69 

1.161626 

-1.173996 

1. 637512 

-65.61 

-1.636781 

6.260937 

6.661765 

lu2 .69 

.90 

-1.692976 

6.662623 

6.652512 

102.97 

1.269613 

-1.270316 

1.  7610  76 

-65.50 

-2.239661 

9. 726229 

9.976769 

102.97 

.99 

-1.663769 

9.666297 

7.„776g3 

193.63 

1.366916 

-  1 . 35o222 

1.911619 

-65. 2C 

-2.671013 

11.166963 

11.501631 

103.63 

.70 

-1.7779J3 

7. 1966c 7 

7.610767 

103.66 

1.633222 

-1.625200 

2. 019091 

-66.90 

-3. 110656 

12.590212 

12.966662 

163.66 

.75 

-1.662616 

7. 371707 

7. bj  63u6 

106.33 

1.669390 

-1.6666-0 

2.091607 

-66.6- 

-3.529902 

13.621950 

16.265571 

106.33 

.00 

-1.961339 

7.361351 

7.612067 

106.76 

1.511925 

-1.675663 

2.112560 

-66. 3C 

-3.662677 

1<*. 720703 

15.226136 

106.76 

,  95 

-1.927322 

7.C7616J 

7.335669 

1C5.23 

1.679933 

-1.627977 

2.055611 

-  6% . 0C 

-6.095558 

15.061053 

15.566676 

105.23 

.90 

-1.796915 

9.386966 

6.632361 

1-5.70 

1.357766 

-1.296632 

1.677593 

-63.66 

-6.037929 

16.366166 

16.922657 

105.78 

.99 

-1.666351 

6.963600 

5.196551 

106.16 

1.075010 

-1.015379 

1.678730 

-63.37 

-3.620636 

11.769025 

12.275309 

106.16 

1.06 

-  •  «iQ3(>3  3 

.  bOCOli) 

.WOOOtl 

105.13 

.000002 

-.63331*2 

.000003 

-61.72 

-.000007 

.000026 

.000 j25 

lw6 . 66 

SECTION  lift  coefficient  a  - 

SECTION  PITCH  HOH  COEFF  a 

SECTION  ROLL  NON  COEFF  ■  - 

X-COOROINATE  OF  SECTION  CP  * 

SPANNISE  COORDINATE  OF  SECTION  CP« 

real 

1.097369 

.963126 

1.769937 

.076665 

2.926692 

INAG 

6.919662 

-.966195 

7.152366 

MAG 

5.060775 

1.368762 

7.363327 

PHASE 

102.57 

-65.63 

103.75 

Program  2  Output 


Unsteady  Loading  Test  Case 
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APPENDIX  F 

PROGRAM  3  OUTPUT  FOR  SAMPLE  CALCULATION 
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Program  3  Output  —  Unsteady  Loading  Test  Case 
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